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NOMENCLATURE 


Drag coefficient based on integrated pressures and skin friction 
(without base drag contribution) 

Drag coefficient including estimated base drag contribution 

Form drag coefficient based on computed momentum thickness at 
"i nfi nity" 

Form drag coefficient based on momentum thickness at the 
trailing edge (Squire-Young formula) 

1 2 

Skin friction coefficient, x/-^ pU 

w 2 e e 

Lift coefficient based on integrated pressure distribution 

Lift coefficient based on circulation determined by "Kutta" 
condition. This is the quantity inputted into the Kutta 
condition in the inverse mode 

Moment coefficient based on integrated pressure distribution 

Pressure coefficient 

Boundary layer shape factor, <5*/e 

Transformed incompressible boundary layer shape factor 
Grid point labels in the w,r directions 
Mach number 

Mach number parameter used in inviscid switch function 
Relaxation factors in the multi -grid/ADI scheme 
Boundary source strength, p g U e 5* 

2 

Parameter appearing in inviscid switch function (Q c = M ) 
Reynolds number 

Arc length along airfoil and wake 

Velocity parallel to body and wake surfaces 

Velocity normal to body and wake surface 

Cartesian coordinates along and normal to chord plane 

Angle of attack 
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NOMENCLATURE (CON'T) 


6 Boundary layer thickness 

5* Boundary layer displacement thickness 

6 Boundary layer momentum thickness 

71 Cole's wake parameter 

^2,3 Parameters controlling artifical viscosity 

p Density 

Skin friction 


SUBSCRIPTS 

B » T Denotes values on lower and upper surfaces of body and wake 

respectively 

C0RR Denotes quantity corrected for wind tunnel interference 

e Denotes boundary layer edge values and inviscid surface values 

te Denotes quantity evaluated at trailing edge 

" Denotes free stream quantity 
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1 . INTRODUCTION 


The computer program, GRUMFOIL, described in this report is for the 
computation of steady, viscous, transonic flow over two dimensional 
airfoils. The code is applicable to flows with freestream Mach numbers less 
than one and high Reynolds numbers for which the boundary layers are thin and 
turbulent over most of the airfoil surface. Our formulation makes use of a 
local analytic solution to correct for strong interaction and normal pressure 
gradient effects at trailing edges. The local solution is based on a formal 
asymptotic analysis for airfoils with cusped trailing edges. Although 
strictly applicable only to airfoils with small included trailing edge angles, 
the present code can be applied to airfoils with arbitrarily large trailing 
edge angles but probably with reduced accuracy of the trailing edge correction 
procedure. The code is intended primarily for use with airfoils with sharp 
trailing edges but can be applied with some approximation to airfoils with 
small trailing edge bluntness. The latest version of the code, which is 
described in this volume, can treat flows with separation and appears to be 
useful up to and including maximum lift. 

The GRUMFOIL code employs a zonal type viscid/inviscid approach to 
viscous flows. It is based on a finite difference solution of the potential 
flow equation for the inviscid flow, integral entrainment methods for the 
laminar and turbulent boundary layers and wakes, and a complete set of 
matching conditions coupling the inviscid and viscid flow fields. A detailed 
description of the basic theoretical and numerical formulation employed in the 
GRUMFOIL code is given in Ref. 1. 

Several versions of the GRUMFOIL code have been developed over the 
years. The preliminary MCMJ-4 version, which was completed in early 1981, 
follows the formulation in Ref. 1. That version employed Jameson's (Refs. 2 
and 3) mixed flow relaxation technique for solving the conservative full 
potential equation using a circle plane ("0" type) coordinate system (A 
nonconservative option was also provided.). The method employed a direct 
solver (Ref. 3) to accelerate convergence. The turbulent boundary layer was 
solved using the standard form of Green's lag-entrainment method. A direct 
iterative scheme was used to treat the coupling between the viscid and 
inviscid solutions. 
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This user manual applies to the most recent version of the code, 
designated MCMJ-9N. This version incorporates several improvements to both 
the numerical procedures used to solve the equations and to the turbulence 
model employed in the lag-entrainment method. 

The improvement to the numerics involved: A switch from the SLOR-Hybrid 

method for solving the potential equations to Jameson's multi -grid ADI method 
(Ref. 4), a more accurate far field treatment, and a more robust second-order 
artificial viscosity relation. The modifications to the lag-entrainment 
method were directed at improving the prediction of separated flows. We also 
replaced the direct method used for the viscous coupling with Carter's semi- 
inverse method, thus allowing for the use of more realistic velocity profile 
relations for separated flow. These changes greatly increased the speed and 
reliability of the code, and significantly reduced the spatial truncation 
error. The new MCMJ-9N version runs about ten times faster than earlier 
versions of the code and achieves the same accuracy with half the number of 
mesh points in each direction. 

The modifications to the turbulent closure relations used in the standard 
lag-entrainment method were directed at improving its predictions of separated 
flows. Integral entrainment methods require the specification of certain 
velocity profile and equilibrium entrainment relations in order to arrive at a 
closed system of equations. In the original Green's lag-entrainment method, 
which was used in the early versions of GRUMFOIL, these were determined from 
empirical correlations for attached turbulent boundary layers. The original 
method was not intended to be used for separated flows, and it is not 
surprising that it did not perform well in the GRUMFOIL code when applied to 
such cases. The present MCMJ-9N version of the code incorporates more 
realistic velocity profile and equilibrium entrainment relations aimed at 
improving the prediction of separated flows. 

The GRUMFOIL code is based on the potential flow equation which ignores 
entropy increases across the shock wave and the associated vorticity 
downstream. Recently, Hafez (Ref. 5) has developed a simple approximate 
procedure for partially accounting for nonisentropic flow effects. In his 
method the entropy effect on the shock relations is accounted for, but the 
effects due to the downstream vorticity are neglected. Results presented in 
Ref. 5 for strong shock cases, where the entropy effect is certainly 
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important, clearly demonstrated the effectiveness of his approach. This type 
of correction is likely to be important in viscous transonic flows near 
maximum lift where strong shock waves and separated flows are present. 
Consequently, we have incorporated the Hafez entropy correction into the 
current version of the GRUMFOIL code. 

All the modifications discussed above, except for the entrainment 
relation changes and entropy corrections, are fully described in Ref. 6. The 
latter improvements will be reported on in a future publication that will deal 
with the extended capabilities of the new code for separated flows. 

In addition to the above revisions, we have also incorporated an improved 
expression for the base pressure which is used in the code to make an estimate 
of the base drag contribution to the total drag. The correction is based on 
Hoerner's correlation for base drag on small rearward facing steps (Ref. 7). 
For completeness, we have included the new expression used for the base 
pressure in the following section. 

The present report is intended to serve as a user's manual for the MCMJ- 
9N version of the GRUMFOIL code. In the following section we provide general 
instructions for running the code and for interpreting the output. In Section 
3, PROGRAM DESCRIPTION, we provide a more detailed description of the program 
which should be helpful for setting up the input data file and for 
understanding how the code functions. A sample calculation including all 
inputs and selected output data is included in the Appendix. 
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2. GENERAL OPERATING INSTRUCTIONS 


In this section we provide some general information that is useful for 
running the MCMJ-9N version of the GRUMFOIL code. We identify all the 
parameters that are needed to set up a computation and provide some 
information that is helpful in interpreting the output data. More detailed 
information that is required to set-up the input data set and to read the 
output is provided in the following section. The present section is divided 
in two parts-one dealing with the input data and the other with the output. 

2.1 Input Data 

Two classes of input data are required - one defines the physical problem 
being addressed and the other establishes the program configuration selected 
for the computation. In most cases the user need only concern himself with 
the first data class which fixes the parameters defining the physical problem 
such as free stream Mach number M^,, Reynolds number Re, and so on. The second 
set of data fixes the program configuration and sets things such as the mesh 
configuration, artificial viscosity, and the velocity profile and entrainment 
relations. Ordinarily the program will be run with these parameters set at 
suggested preset values which will be appropriate for most problems except 
near maximum lift where there are extensive regions of separation present. 

2.1.1 Problem Definition 

This information defines the problem being addressed. The user must 
first specify the value of the five constants: 

IBL, LMODE, Ma,, Re and either a or C L . 

The parameter IBL determines whether the computation is to be inviscid (IBL=0) 
or viscous (IBL=1) and the constant, LMODE, determines whether it is to be a 
direct, a prescribed, (LM0DE=0) or inverse, C L prescribed, (LM0DE=-1) 
computation. The user is cautioned not to run in the inverse mode at or near 
maximum lift. If a is prescribed that is greater than maximum lift no 
solution exists and the code will fail. The quantities M^ Re, a, and C L are 
the free stream Mach number, chord Reynolds number, angle of attack and lift 
coefficient respectively. The problem definition is completed by specifying 
the airfoil ordinates, Y(X), 
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which need not be normalized by any particular reference length and by giving 
the value of the trailing edge included angle (in degrees), 0^, and the slope 
of the trailing edge bisector. When neither is known accurately, their true 
values will be computed within the program from the tabulated ordinates if 
Sj-g is given a dummy negative value. 

The program has been written to permit a sequence of problems to be 
solved for a specified set of values of M ro and a or C^. To exercise this 
option the user must specify the ranges (M^, M 2 ), (C^, C^) or ^ a1 ’ a2 ^ 017 
parameters over which solutions are sought and the increments zC|_, Aa and 
with which the program is to increment the values of M„, a or C|_. 

The parameters identified above provide the minimum data that needs to be 
specified to carry out a given calculation. Although the remaining parameters 
to be discussed below can usually be left at their preassigned values there 
may be occasions where it is desirable to change them. The following 
discussion is intended as a guide to help understand the rationale to be used 
in deciding when and how these values should be changed. The remaining 
parameters essentially control the numerical and theoretical formulation to be 
used in the computation. The program configuration parameters fall into four 
main groups: 1) inviscid solution, 2) boundary layer solution, 3) 
viscid/inviscid relaxation parameters, and 4) convergence tolerances. 

2.1.2 Inviscid Solution 

The parameters in this group control the numerical method used to solve 
the inviscid potential flow equation. The Hafez entropy correction is 
controlled by the parameter IRH and is either 0 or 1, 

irm [ 0 No entropy correction 

‘ 1 Entropy correction applied 

The default value is IRH=0 so that the entropy correction is not applied in 
the baseline version. The meshes employed in the inviscid calculation are 
controlled by the parameters, 

FNC, FNR, FMESH and 6MESH 

The potential flow calculation is ordinarily carried out on a sequence of 
macro-meshes. A complete solution is first carried out using multi-grid, on 
the coarsest macro-mesh. The parameters FNR and FNC specify the number of 
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mesh divisions in the R and w directions on the coarse grid respectively, 
where R, ware polar coordinates in the circle plane. The multi-grid method 
employs a sequence of coarser "micro" meshes that are formed by deleting grid 
points from the preceding micro-mesh. The parameter 6MESH, specifies the 
number of micro-grids employed in the multi -grid solution on each macro-grid, 
and must be specified for each macro-grid. The value of GMESH must be chosen 
so that the coarsest grid used in the multi -grid solver contains at least two 
divisions in R. The parameter, FMESH, specifies the number of macro-grids to 
be employed. If FMESH is greater than one, the solutions carried out on the 
coarser grids will be used to generate initial data on the succeeding macro - 
meshes. In the baseline configuration the code is set-up to use three macro- 
grids (FMESH=3), starting from a (40x8) Mesh (FNC=40 and FNR=8) . Thus, the 
computations will be carried out on a sequence of three grids (40x8), (80x16) 
and (160x32) with the multi -grid parameter GMESH set to 3., 4., 5., respec- 
tively. In order to avoid runaway divergence in an inviscid calculation a 
maximum number of multi-grid cycles must be assigned. Typically, 25 is used 
for the parameter T0T0 on each macro-mesh. In a viscid calculation the 
boundary layer is updated after each complete multi -grid cycle, therefore, 

T0T0 is always set to 2. 

The artificial viscosity is controlled by three parameters VIS1, VIS2, 
and VIS3 which correspond to the parameters \^, \>>, vj used in Ref. 6. The 
quantity VIS2 determines the order of the scheme - second order if VIS2=0 and 
first-order if VIS2=1. The quantities VIS1 and VIS3 set the levels of the 
first- and second-order artificial viscosities respectively. The baseline 
program is set to run second-order accurate (VIS2=0) with VI Sl=l and V IS3=2. 
The switch function which turns off the artificial viscosity at subsonic 
points is controlled by the parameter Q c . The value set for Q c determines the 
local Mach number, M c , (Q C =M C ^) at which the artificial viscosity is reduced 
to zero. We recommend a default value of Q c = 0.90. 

The Parameters, (i=l 5) are relaxation parameters employed in the 
ADI scheme. The constants Pj, P2, P3 are proportional to the 0^, c^, 
constants in Ref. 4 that set the so-called c^, ^ and ^ temporal damping 
terms in the iterative operator. The P4 and P5 terms are equal to the 
relaxation factor multiplying the residual at supersonic and subsonic points 
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respectively. The relaxation values will usually not have to be changed from 
the default values set in the baseline code, which are 

Pi = .25 P 4 = .60, .70, .80 

P 2 = .45 P 5 = 1.5 

P 3 = .45 

The three values of P 4 are the values on the three macro grids (40x8), (80x16) 
and (160x32). The other P^ can generally be kept the same on the three macro- 
grids. 

The last parameter that must be set in the inviscid calculations is NINV 
which is associated with an iteration that determines the angle of attack from 
the Kutta condition when the program is operated in the inverse mode (LM0DE = 
-1). In this mode the Kutta condition reduces to a non-linear algebraic 
equation which must be solved by iteration to determine the angle of attack. 
The parameter NINV sets the number of cycles that are to be carried out in 
this computation. The default value is NINV = 4. This is sufficient for 
calculations carried out to maximum residual levels of 10"® to 10"^. The 
value of NINV should be increased if higher levels of convergence are sought 
in the inviscid computations. 

2.1.3 Boundary Layer Solution Procedure 

This group of parameters controls the formulation used for the boundary 
layer computation. The quantities ICURV and IDIV control the terms appearing 
in Green's lag-entrainment method for the effects of curvature and mean 
dilation on turbulence. These terms are included if ICURV and IDIV are set to 
one (default) and are dropped if they are set to zero. The present code 
provides some options on the choice of the velocity profile, equilibrium 
entrainment and skin friction relations that can be used with the lag- 
entrainment method. The velocity profile relation is controlled by the value 
of the parameter, IH1HB. The various possibilities are described in Ref. 6. 


IH1HB=1 

relation in original Green's 
lag-entrainment method 

IH1HB=2 

East, Smith 

and Merryman 

IH1HB=3 

Le Balleur 
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The choice of the equilibrium entrainment relation is controlled by the 
value of the parameter IMCE. Two options are provided. The value IMCE=1 
selects the equilibrium entrainment relation employed in the original Green's 
lag-entrainment relation. As mentioned previously this relation greatly 
overestimates the levels of entrainment in separation zones. The choice IMCE 
= 2 employs an equilibrium entrainment relation proposed by Le Balleur (Ref. 

8) for separated flows. It is based on a smooth blend between Michel's (Ref. 

9) entrainment relation for attached flow and a mixing layer type formula for 
separated flow. 

Although these various formulations for the velocity and equilibrium 
relations are clearly improvements to the original method when the flow is 
separated we do not yet know which is the best one to use in general. On the 
basis of limited calculations we recommend, at this time, the values IH1HB=3 
and IMCE=2 which have been set as the default values in the code. 

The skin friction relation employed in the original Green's lag- 
entrainment method does not produce realistic predictions of the skin-friction 
levels in regions of reversed flow. Since the contributions of the skin 
friction to the momentum balance is not important in separated flow zones, the 
skin friction was set to zero in such regions in the original GRUMFOIL code. 

An improved skin friction correlation that produces reasonable predictions of 
skin friction in regions of ‘ reversed flow was proposed by Swafford in Ref. 10 
on the basis of limited separated flow data. The present version of the code 
allows for the choice of using either the original Green's or Swafford's skin 
friction relation in the turbulent boundary layer solution. This option is 
controlled by the value of a parameter ICF. The original Green's skin 
friction relation will be employed if ICF=1 while Swafford's will be used if 
ICF=2. A default value, ICF=2, using Swafford's skin friction relation has 
been set in the baseline version of the code. 

The user must also specify how the transition point locations are to be 
determined. These are controlled by the values of the parameters XTRANU, 
XTRANL, and KT where the quantities XTRANU and XTRANL are the transition point 
locations in percent chord on the upper and lower surfaces of the airfoil, 
respectively, and KT is a parameter that specifies the method for determining 
the natural transition point locations. The choices KT=1, 2, or 3 will select 
the natural transition criteria of Crabtree (Ref. 11); Michel (Ref. 12); or 
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Stevens, Goradia, and Braden (Ref. 13); respectively. In these cases the 
parameters XTRANU and XTRANL should both be set to zero. If the user wishes 
to directly assign the location of the transition points on the upper and 
lower surface of the airfoil, he should set KT=4 and set the parameters XTRANU 
and XTRANL to the desired locations in percent chord. Finally, if KT=1, 2, or 
3 and XTRANU and XTRANL are assigned non-zero values, the program will select 
the most forward of the two possible transition positions. Since KT=1 and 
XTRANU, L = .10 are defaults, this is the default mode of code operation. 

The program will also check for laminar separation and will use 
Crabtree's criteria (Ref. 11) to decide if it is a short or long bubble 
leading-edge separation. In either case, all other transition criteria will 
be bypassed and the transition point will be set to the laminar separation 
point position. This procedure is physically sound only if the leading edge 
separation is of the short bubble type. For long bubble separation, the 
procedure, though not physically correct, will assure that a solution will be 
completed. 

We have found it desirable to override the user specified input values of 
XTRANU and XTRANL under two circumstances. First, if the selected transition 
point location is upstream of the stagnation point the code will automatically 
reset the assigned transition point to the trailing edge and the calculation 
will continue effectively as if in the natural transition mode. This is a 
physically sensible procedure since transition to turbulent flow should not be 
possible if a transition point location is chosen to be upstream of the 
stagnation point in a strong, favorable pressure gradient. The second set of 
circumstances where we have found it desirable to alter the assigned 
transition point locations is where, although the transition point is 
downstream of the stagnation point, it is too close to the stagnation point. 

We have found that the code will frequently fail in the turbulent boundary 
layer solution if transition is assigned too close to the stagnation point in 
a strong, favorable pressure gradient. Therefore, a Reynolds number based 
upon the momentum thickness at the assigned transition point is calculated; if 
it is less than 225, we reassign the transition point to the first grid point 
downstream where this Reynolds number is greater than 225. In both of the 
cases described above, a notification message is printed in the tabular 
output. 
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In order to prevent run-away divergent solutions we have set a maximum on 
the number of global viscid-inviscid cycles. This is controlled by the 
parameters, ITBLMO, ITBLM1, and ITBLM2, which specify the maximum number of 
global cycles on each mesh. These have been given default values 22, 22 and 
30 on the coarse, medium and fine meshes, respectively. They can be increased 
if further cycles are needed in slowly convergent solutions as is usually 
required when extensive separation is present. 

2.1.4 Viscid-inviscid Relaxation Parameters 

The global viscid-inviscid calculation requires significant 
underrelaxation to achieve convergence. Separate relaxation factors are 
introduced into the solutions for the source strength on the airfoil, the 
source strength in the wake and for the wake curvature terms. These are given 
by RQE, RQW and RBW, respectively. The values of these relaxation factors can 
be set to different values on each of the macro-grids. The interaction on 
each macro-grid starts with the values assigned to RQE, RQW and RBW and are 
slowly increased to larger values as the iteration proceeds. The final values 
are multiples of the initial values given by the factors WRQ, WRQW and WRBW. 

We must also set values of an auxiliary relaxation factor RTE which controls 
the iteration of a constant that appears in the wake thickness coupling 
conditions. The default values set for these parameters are given in the 
following table. 

COARSE GRID (40X8) 


WRQO =5.0 

RQEO = 0.2 

WRQW0= 5.0 

RQWO =0.2 

WRBW0= 5.0 

RBWO = 0.2 


RTEO = 1.0 

MEDIUM GRID (80X16) 

WRQ1 = 1.0 

RQE1 = 1.0 

WRQW1= 1.0 

RQW1 = 1.0 

WRBW1= 1.0 

RBW1 = 0.4 


RTE1 = 1.0 
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FINE GRID (160X32) 



WRQ2 = 1.0 

RQE2 = 0.8 


WRQW2= 1.0 

RQW2 = 0.4 


WRBW2=1 .0 

RBW2 = 0.20 
RTE2 = 0.10 

These 

values produce rapid and sure 

convergence for attached flows with shock 

waves 

that are not too strong. For 

more difficult problems with extensive 

separation, use the following reduced values (shown underlined): 


COARSE 

GRID (40X8) 


WRQO = 

RQEO = 0.2 


WRQW0= 2^5_ 

RQWO =0.2 


WRBW0= US_ 

RBWO = 0.2 
RTEO = 1.0 


MEDIUM GRID (80X16) 


WRQ1 =1.0 

RQE1 = 0.5 


WRQW1= 1.0 

RQW1 = 0.4 


WRBW1= 1.0 

RBW1 =0.2 
RTE1 = 1.0 


FINE GRID (160X32) 


WRQ2 = 1 .0 

RQE2 = 0^ 


WRQW2= 1.0 

RQW2 = 0^ 


WRBW2= 1.0 

RBW2 = 

RTE2 = 0.05 


The use of these reduced values will help convergence but at the expense of 
longer run times. 

2.1.5 Convergence Tolerances 

The solution is terminated when certain monitored quantities are reduced 
below prescribed tolerances. The four quantities monitored during the 
calculations are: 1) the maximum residual, 2) the change in angle of attack 

or lift coefficient, 3) the change in trailing edge pressure, and 4) the 
maximum difference between the boundary layer and inviscid solutions for the 
inviscid streamwise velocity on the airfoil surface and in the wake. The 
parameters controlling the tolerances with the recommended default values are 
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COV = 0.5 X 10 " 4 
COVA, COVL = 0.005 

COVP = 0.001 
DUCR = 0.01 

Terminating the solution 
the angle of attack, 
places, 


Maximum residual 

Change in angle or attack or 
lift coefficient 

Change in trailing edge pressure 

Difference in viscous and inviscid 
streamwise velocity 


at these values will usually produce solutions with 
lift coefficient and pressures converged to three decimal 
and the drag coefficient to four places. 


2 . 1.6 Experimental Data 


The program allows the user to input experimental data for the pressure 
distribution Cp(x); displacement and momentum thickness <5*(x) and 0 (x); and 
skin friction coefficient, Cf(x). These data will be incorporated into the 
final solution output plots. 

2.1.7 Summary Recommendations For Difficult Cases 


The program has been organized so that it will run with a minimum amount 
of user supplied data. Ordinarily, for most cases the user will have to 
supply only data that define the physical problem being addressed. The data 
which control the program configuration has been set at default values that 
provide a reasonable compromise between convergence speed and reliability. 

The default values should be satisfactory for all cases that result in 
attached or mildly separated flow. Adjustment of some of the program 
configuration parameters will be necessary for cases near maximum lift which 
involve extensive regions of separated flow. In the following, we summarize 
the changes that we have found to be effective for these more difficult cases. 

If, for supercritical cases, the inviscid solution shows osciliatory 
convergence or divergence as evidenced by the plot of convergence history 
included in the output, this is usually an indication that the ADI relaxation 
factor, P 4 , is too large. Convergence can be improved by reducing the value 
of P 4 by a decrement of 0.1 or 0.2 on the offending grid. 

If there is extensive separation, it is always necessary to use the 
smaller set of global relaxation factors given earlier. In these cases the 
global viscid-inviscid interaction will be significantly slowed thus requiring 
the program to run longer to achieve the same level of convergence. There- 
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fore, the parameters ITBLMO, ITBLM1, and ITBLM2, which set the maximum number 
of global iterations on each of the three macro-meshes, should be increased to 
100, 150 and 150, respectively, to allow for the larger number of iterations 
requi red. 


We have also found that, in the separated flow cases near maximum lift it 
is very helpful to start the calculation on the next finer macro-mesh (80 x 
16). Under these circumstances only two macro-meshes are used (FNC = 80, FNR 
= 16, and FMESH = 2) and it is necessary to use the reduced global relaxation 
factors appropriate to the finer mesh start as follows: 


WRQO =2.0 
WRQWO =1.5 
WRBWO = 1.0 


WRQ1 = 1.0 
WRQW1 = 1.0 
WRBW1 = 1.0 


FIRST GRID (80X16) 
RQEO = 0.15 
RQWO = 0.10 
RBWO =0.10 
RTEO = 1.0 

SECOND GRID (160X32) 
RQE1 = 0.2 
RQW1 = 0.2 
RBW1 = 0.1 
RTE1 = 0.05 


These recommendations are intended only as helpful guidelines for 
treating separated flow cases near maximum lift. The user may have to 
experiment with the setting of these parameters in order to obtain successful 
solutions for these cases. 


2.2 Output Data 

Various options are provided to the user in selecting the amount of 
output data desired. The procedure for controlling the output is described in 
the following section. The principle mode for data output is in the form of 
graphs of the pressure distribution, the distribution of various boundary 
layer quantities, and a plot of convergence history. For completeness, the 
values of the forces, moments, and other numerical data of interest are also 
included in the plots. 


The graphical output data consists of six plots. (A seventh plot is 
produced if experimental pressure data is to be plotted.) The first plot 
contains two surface pressure distributions; one is the "outer" inviscid 
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solution and the other is the viscous composite solution. The difference in 
these two pressure distributions is indicative of the pressure correction 
across the boundary layer. The numerical data on the bottom of the plot 
includes the solution for the lift, drag, and moment coefficients, C|_, C D , and 
C|v| as computed from an integration of the composite pressure distribution. 

Also included are: a title identifying the airfoil geometry; the free stream 

Mach number, M; the angle of attack, ALF; the "circulatory" lift coefficient 
CLC; the number of inviscid iterations carried out on the particular grid 
corresponding to the plot, NIT; the number of mesh division (NC x NR) employed 
in the w; and R directions, the value of the maximum residual in the inviscid 
solution at convergence, RSID; the value of the trailing edge pressure, CPTE, 
and a drag coefficient, CDB, which includes a contribution for the base drag 
of blunt based airfoils. The circulatory lift coefficient, CLC, is the value 
of lift computed from the Kutta condition in the direct mode ( a prescribed) or 
inputted into the Kutta condition, (C|_ prescribed). The value of CLC should 
agree with the value determined from the integration of the outer inviscid 
pressure distribution around the airfoil in limit of zero mesh. We have found 
agreement to within 0.005 on the 160 x 32 grid. 

The base drag contribution to the total drag is obtained from an 
adaptation of Hoerner's correlation (Ref. 7) for rearward facing steps. In 
using Hoerner's result to compute the base drag on airfoils with a small 
trailing edge thickness we have assumed that one of the dividing streamlines 
could be replaced by a solid surface. For completeness we write out the 
expressions used to compute the total drag including our approximation for the 
base drag. It can be written in the form. 


CDB = CD - £ d 

where the base drag contribution ^Cg is given by 

9 9 3.5 

M t c 1 + 0.2 

S = t C pte “ 73T~ 7~~rrr. 2 Cp B ^te cos a 


M 


1 + 0.2 M 


te 


where Cp te , and M te are the trailing edge pressure and Mach number as computed 
from the viscous composite solution, H te is the base height of the airfoil in 
percent chord, a is the angle of attack, and Cpg is the incompressible base 
pressure coefficient normalized by trailing edge values from the viscous 
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composite solution. The base pressure coefficient, Cp B , is determined from 
Hoerner's correlation. It is written in the form: 

Cp B = 0.231 {H te f&) l/3 

where 6 is a boundary layer thickness parameter defined by 

6 = W + ^in = ^TOTAL 

The quantities <^ ax and are the larger and smaller values of the boundary 
layer thickness at the trailing edge on the upper and lower sides of the 
airfoil, respectively. These results are consistent with Hoerner's 
correlation in the limits ^-jn/^ax = 0 and ^ip/^ax = * within the 
approximation that the dividing streamline can be replaced by a solid 
surface. The above expression provides a simple linear interpolation in 
Vin^max between the two limits. The expression for the base drag is 
intended as a simple and very rough estimate, and is considered an improvement 
over the expression used in the earlier versions of GRUMFOIL which neglected 
the contribution of the base pressure term, Cp B . 

The second plot provides a graphical presentation of the iteration 
history. The figure includes plots of the trailing edge pressure and the 
logarithm of the "error" in the inviscid solution. The "error" is defined as 
the ratio of the average residual to its initial value on the particular grid 
for which the results apply. The unit of "work" is defined to be equal to the 
equivalent number of SLOR steps. Information is included at the bottom which 
gives the value of the average residual at the start and finish of the 
computation, RESID1 and RESI02 respectively; the number of inviscid cycles, 
NIT, the total value of work at the end of the computation; and the overall 
convergence rate which is defined by 


RATE 


/RES1D2 i 

tenrr J 


1/Work 


The third plot presents the solution for the displacement thickness, 6 , 
and shape factor, H, on the upper and lower surfaces of the airfoil and in the 
wake. Numerical data at the bottom of the plot include the values of free 
stream Mach number; angle of attack; the transition option employed through 
the value of KT; and the actual transition point locations employed in the 
computations on the upper and lower surfaces of the airfoil, XTRANU and 
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XTRANL, respectively. Also included are the values of the skin friction drag 
and form drag. The quantities CDF, CDFT and CDFB are total skin friction drag 
and the contributions from the upper and lower surface boundary layers, 
respectively, as computed from integrations of the skin friction distribution. 
Two values of the form drag are given: one, COSY, computed from the Squire- 

Young formula, and the other, CDINF, determined from the computed momentum 
thickness at the last point in the far wake. The contributions from the 

momentum thickness of the upper, CDIT, and lower, CDIB, sides of the wake are 

also given. The form drag should equal the total drag in subcritical 

(shockless) flow, and should equal the total drag less the wave drag in 

supercritical flow with shock waves. 

The fourth plot presents the solution for the momentum thicknesses and 
skin friction distribution on the upper and lower surfaces of the airfoil and 
wake (The skin friction is assumed to be zero in the wake.). As discussed 
previously, the skin friction is set to zero in separated flow regions if 
ICF=1 . 

The fifth plot presents the solutions for the entrainment coefficient, 

CE, and the equivalent source velocity which is used as a boundary condition 
in the inviscid outer solution. On the airfoil surface, 0 < x <1, the 
quantity, V, represents the equivalent source velocity on the upper and lower 
sides of the airfoil while in the wake, x > 1.0, it represents the average of 
the source velocity on each side of the wake streamline. The sixth plot 
repeats the viscous solution for pressure distribution and numerical data 
given in the first plot. 

If experimental data for the pressure distribution is included in the 
input data, a seventh plot will be generated that compares the experimental 
pressures with the viscous composite solution. As described in the next 
section, the program is structured to make it convenient to specify the plot 
scales desired. We also permit the user to input a blockage correction in the 
form of a corrected free stream Mach number. If the program is run at a 
corrected Mach number, and the experimental Mach number is included with the 
input of experimental data, the program will calculate blockage corrections to 
the experimental pressure distribution and lift and drag according to the 
formulas 
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3.5 


When running the code in the inverse mode with the lift coefficient 
chosen to match the experimental lift coefficient, it is necessary to make a 
slight adjustment to the input lift coefficient. This small correction is 
necessary because the lift coefficient computed from the integrated pressure 
will differ slightly from the circulatory lift coefficient, CLC, which is the 
quantity input into the Kutta condition in the inverse mode. If we wish the 
theoretical integrated lift coefficient to agree with the experimental value 
of the lift coefficient, we must slightly adjust the value of the input value 
of CLC. The appropriate increment is usually very small and can be quickly 
estimated by trial and error. The increment DCLC, if included in the input 
data, is written at the bottom of this plot. Its default value is 0.0. 


The graphical data is available in tabular form along with tables of 
other properties of the solution that may be of interest. The procedures for 
selecting the desired output options are described in the following section 
along with further details for interpreting the tabular output data and for 
setting up the input data deck. 
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3. PROGRAM DESCRIPTION 


3.1 Program Structure 

This subsection contains a note on compiler and linkage information 
followed by a list of subroutines (with a brief description of each), a 
subroutine tree diagram and a simplified flow chart of the program logic. The 
program is written in FORTRAN 66 and requires each do loop to be executed at 
least once. 

3.1.1 Compiler & Linkage Information 

At Grumman, the program is run in double precision on the IBM 370/3081 by 
employing the AUTODBL feature of the Fortran H extended compiler. This 
recommended precision is about the normal level for CDC machines. The program 
requires about 1000 K bytes of storage in the mode run at Grumman. Some 
installations may require overlay. 

3.1.2 Subroutine Glossary 

We now list the principle subroutines used in the program along with a 
very brief definition of their function. In subroutines GRAPH, CVPLT, and 
AUXPLT are calls to a number of Cal comp routines which are members of a 
Grumman library. The user will probably have to modify these routines 
slightly to make them compatible with the user's graphics installation. The 
normal operation of subroutine RELAX produces a large number of underflows. 

We use subroutine ERRSET, a member of our IBM Fortran library, to suppress 
their write out. The user should find an equivalent replacement. 

COMAP Transforms airfoil into unit circle. 

SPLIF Cubic spline fit routine. 

C0NT0 Determines arc length, slopes, and curvature of the 

profi le. 

ARCL Calculates arc length using the first three derivatives of 

x and y. 

PARAF Parabolic fit routine. 

INTPL Interpolation using Taylor Series. 

SURFL Constructs the derivative of arc length by conjugation. 


18 



CONJ 

FOUCF 

FFORM 

BLDATA 

CPXC 

DISC 

FARFL 

ESTIM 

AIRFL 

LE 

INVBL 

CUT 

DGAM 

DELTAJ 

RELAX 

STAG 

XPAND 

ADD 

VELO 


Conjugation by fast Fourier transform. 

Determines the Fourier coefficients by fast Fourier 
transform. 

Fast Fourier transform. 

Reads and writes the values of the parameters for the 
viscous calculations. 

Corrects experimental values of Ci , Cp, and C„ to match a 
change in indicated experimental Mach number, and scales 
variables for plotting. 

Calculates the mapping modulus exterior to the unit circle 

Calculates the far field functions. 

Initializes the reduced velocity potential to zero. 

Evalautes airfoil coordinates at circle plane nodal points 
(equal circle angle increments). 

Calculates circle plane index of leading edge. 

Provides an initial parabolic distribution of source 
velocity over the rear portion of the profile (on both 
surfaces) to precondition the first inviscid solution. 

Calculates the coordinates of the cut (6= 0° in circle 
plane) and its local slope. 

Computes the change in circulation along the wake due to 
its curvature. 

Calculates function used in determining jump in potential 
across wake more accurately. 

Solution of the flow equations in divergence form using 
rotated difference scheme - Relaxation step. 

Uses the stagnation point location to determine the index 
for the start of the upper or lower surface relaxation 
sweep on every mesh. 

Reorders certain macro-mesh arrays to conform to each 
micro-mesh expansion or contraction. 

Adds correction to the solution on a finer mesh using 
second-order interpolation and smooths the result. 

Calculates surface velocities, Mach numbers, and pressures 



QWK 


Calculates the velocity on each side of the wake and the 
ordinates of the wake streamline. 


WAKE Forms basic arrays for edge velocity and coordinates of the 

profile and wake for boundary layer calculation. 

JUMPC Calculates a surface shock position as the last upstream 

point in a supersonic region. 

JUMP Calculates the pressure downstream of a shock treated as an 

isentropic compression and also from the Rankine-Hugoniot 
relations. 

XTRAP Determines trailing-edge values of functions by linear 

extrapol ation. 

FORCS Calculates force and moment coefficients from Simpson's 

Rule integration of pressures in circle plane. 

CPLOT Generates a crude pressure plot around airfoil and into 

wake. 

CHART Generates a Mach number distribution chart for the whole 

flow field. (Can be made to generate similar charts for 
the total velocity or either of the component velocities by 
changing the last parameter in the calling statement.) 

SONICA Stores a switch function along each radius of the 

computational plane for use in the calculation of the sonic 
line. 

SONICB Calculate the coordinates of the sonic line. 

SONIC 

GRAPH Cal comp plotting package for C p , 5*, H, e, Cf, C^, and V a 

CVPLT (source velocity) plots. Also^for the convergence history 

AUXPLT of the average residual and trailing edge pressure. Also 

for auxilliary C p plot of assignable scale to be used in 
comparisons with H experiment. 

VISCAL Calculates the boundary layer over the airfoil and into the 

wake. Uses a direct solution for the laminar run and an 
inverse solution for the turbulent flow. Corrected source 
velocities are computed from the updated source strength 
function and introduced into the boundary conditions of the 
inviscid solution. 

VISID Initializes the laminar and turbulent boundary layer 

calculations. 

THWAIT Calculates laminar boundary layer by Thwaite's integral 

method, tests for transition, and generates initial values 
of Cf, H, and 6* for the turbulent boundary layer 
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calculations. Also integrates Cf to get the laminar 
contribution to the skin friction drag. 

BLKDATA Initializes arrays of Thwaite's functions for laminar 

boundary layer calculation. 

EQENTO Calculates equilibrium entrainment coefficient. 

EDGEC Uses linear fits to the boundary layer source strength/edge 

velocity and curvature to calculate their derivatives which 
are used in interpolation. 

RKGM Runge-Kutta-Gil solution of Green's modified lag- 

entrainment turbulent boundary layer equations. 

GREENC Evaluates the right hand side of Green's equations which 

are written in terms of circle plane coordinates and source 
velocity. 

OUTBOC Loads the arrays for UEV, Cf, S*, and eat circle plane 

nodal points. 

STOX Prints out table of boundary layer characteristics as well 

as UEI and UEV. It also prints out the maximum difference 
between UEI and UEV and its position as well as the average 
difference. 

PIWIG Iterative solution for the boundary layer thickness from 

the skin friction and displacement thickness using law of 
the wal 1 . 

G12 Solutions of the universal trailing edge correction 

functions. 

FCAP Calculates the flow deflection angles along the wake. 

SMTH Uses running average smoothing for uniformly spaced 

function. 

BFIT Assigns parabolic variation of the flow wake angle near the 

trailing edge with zero at the trailing edge. 

SPREAD Linear interpolation of expanded arrays at mesh division. 

3.1.3 Subroutine Calling Sequence 

The accompanying Subroutine Tree (Fig. 1) shows how the various program 
subroutines are called. In addition, the subroutines which are related 
specifically to the viscous calculation are enclosed in rectangles. 
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Fig. 1 Subroutine Tree 
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3.1.4 Program Flow Chart 


The simplified flow chart in Fig. 2 gives a general overall view of the 
program. The main features to be noted are: 

o The IBL branches where we switch to either a viscous calculation (IBL 
= 1) or an inviscid calculation (IBL = 0). 

o The ADI branch for expanding and contracting meshes. 

o The convergence check on the relaxation procedure. 

o The looping back for boundary layer iteration. 

o The looping back for mesh division. 

o The looping back after incrementing Mach number, angle of attack or 

lift coefficient. 

3.2 Input Data File 

This section lists and defines all essential parameters for executing 
program GRUMFOIL (MCMJ-9N). They are given in the order in which they are 
used in setting up a calculation. It should be emphasized that the majority 
of listed parameters in the program NAMELISTS rarely need to be changed from 
their default values. At the end of this section. Fig. 3 depicts the makeup 
of an input data file as an equivalent card image file. The data file used to 
generate the sample case is listed in the Appendix. 

A number of header cards are placed in the data file to help in its 
organization and as cues for the user. They are necessary parts of the data 
file and must be included. The labels shown on them, identifying the values 
of the parameters which appear on the following card or cards, match the 
variable names read by the program. These names are defined below and the 
pertinent program read and format statements are also included. 

The first entry in the file is used to describe the run. A maximum of 80 
characters can be used for an alpha-numeric description of the profile and 
flow situation being analyzed. The statements in the program are 

READ( IREAD,530) (TITLE(I) ,1=1,20) 

530 F0RMAT(20A4) 

In the above, IREAD represents the system input data set reference 
number. Similarly, within the code IWRIT is used to designate the output 
number. These parameters are set at the beginning of the main program and 
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Fig. 2 Simplified Flow Chart 
































must be modified to comply with the system being used. The next record is 
used as a header card to identify the items to be read from the following 
record. This first dummy entry as indicated from the sample case shows FNC 
FNR and FMESH. It is read from the program by: 

RE AD (I RE AD, 500) 

500 FORMAT (IX) 

The next record contains the actual values of these three parameters and is 
read from the program by: 

READ(IREAD,510)FNC, FNR, FMESH 
510 FORMAT (8F10.6) 

where: 

FNC: REAL Cols 1-10 Number of angular divisions on the first, 

coarse macro-mesh. Suggested minimum 40., 
maximum on the final macro-mesh is 192. 

FNR: REAL Cols 11-20 Number of radial divisions on the first, 

coarse macro-mesh. Suggested minimum 8., 
maximum on the final macro-mesh is 32. 

FMESH: REAL Cols 21-30 Number of macro-meshes to complete 

calculation. Usually 3. 

The next record is used to establish the program mode. The pertinent 
statements in the code are: 

READ ( IREAD ,661) ICASE , IBL ,LM0DE , IOPUT 
661 FORMAT (11,14,515) 

where: 

ICASE Integer 1 -*■ 9 Col 1 Used to identify run or data set 

number. 

IBL Integer 0 or 1 in Col 5 0 inviscid mode 

1 boundary layer included 

LMODE Integer 0 or -1 in Cols 9-10 0 direct mode; a prescribed, find C, 

-1 inverse mode; C^ prescribed, find a 

IOPUT Integer 0 -*■ 7 in Col 15 IOPUT is one of two parameters 

used to control the printed and 
graphical output of the code. 

This one effects material which is 
produced for both inviscid and 
viscous calculations. It is also 
used to print some material which 
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TABLE I - Output Generated for Various Values of IOPUT 


IOPUT -*• 
OUTPUT + 


Final Mesh Cal comp Plots. 

♦Convergence Summary Table Each Mesh. 


Airfoil Geometry. Final Mesh: convergence, 

shock position, forces, surface pressure 
table, Richardson extrapolation of forces. 


All Meshes: convergence, shock 

position, forces, 

♦wake pressures and circulation for last 
boundary layer iterate. 


Mapping Convergence All Rough C_ Plots. 
♦Convergence, shock position, forces for 
every boundary layer iterate. 


♦All Meshes: body and wake curvature for 
last boundary layer iterate. 


Mapping Coefficients. Intermediate Mesh 
Cal comp Plots. 


Final Mesh Mach Number Chart. 


Intermediate Meshes: surface pressure 

table, Mach number chart, 

♦for every boundary layer iterate. 








only applies to a viscous run. 

The accompanying Table I indicates 
the increasing amount of material 
generated by increasing values of 
IOPUT. The items marked by an 
asterisk only apply to viscous 
cases. We recommend the value 1 
for this parameter. The second 
parameter, IOUT, which controls 
boundary layer output is described 
1 ater. 

The next record or continuations thereof refers to a NAMELIST statement 
with 25 arguments. The record must be included but need not have any entries, 
in which case, all default values will be used. The program read statement 
is: 


READ( IREAD,MXF) 


where the parameters are: 


Parameter 

NINV 


IRH 

COVA 


COVL 


COVP 


Comments Default 


Integer, number of iterations used to determine 4 
angle of attack from the Kutta condition when 
program is operated in the inverse mode (LMODE = -1). 

Can be doubled if convergence levels in an inviscid 
case below 10" 1U is sought. 

Integer, Hafez entropy correction, = 0 no entropy 0 
correction, = 1 entropy correction applied. 

Real, convergence criterion applied to angle of .005 

attack when program is run in the inverse mode. 

If this criterion and that for the relaxation 
convergence, and COVP, and DUCR are met, no 


further boundary layer iterations are done on the 
particular mesh. 

Real, convergence criterion applied to lift .005 

coefficient when program is run in normal 

mode. If this criterion and that for the 

relaxation convergence, and COVP, and DUCR are met, 

no further boundary layer iterations are done on 

the particular mesh. 

Real, convergence criterion applied to the .001 

trailing edge pressure. If COV, COVA or COVL 
and DUCR are also satisfied, no further boundary 
layer iterations are done on the particular macro- 
mesh. 
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Parameter 


Comments 


Default 


OUCR Real, convergence criterion applied to the .01 

maximum difference between UEI, the in viscid 
velocity distribution, and UEV the edge velocity 
from the boundary layer equations. Check is made 
on upper and lower surfaces and in the wake. 

If COV, COVA, or COVL and COVP are also 
satisfied, no further boundary layer iterations 
are done on the particular macro-mesh. 


ITBLMO, Integer, maximum number of boundary layer 

ITBLM1, iterations on successive macro-meshs. They are 

ITBLM2 meant to limit time spent on poorly converging 

cases. In cases with large separation zones, 
when reduced relaxation factors are called for, 
increasing these parameters (e.g., 100,150,150) 
may be required for satisfactory convergence. 

The remaining 16 NAMELIST parameters are designed to be used in 
comparisons with experimental data. Six are used to generate an auxiliary 
scaled pressure plot of results generated by the code. The user selects the 
scale as detailed below (defaults produce the same scale as used by the 
code). Seven are used when experimental data are available to add to the 
auxiliary pressure plot. A blockage correction will automatically be applied 
to these data when the calculation Mach number is different from the 
experimental Mach number. The last two are used when experimental boundary 
layer data are available to add to the 6*, 0, H, and C f plots. 

IPSW Integer, = 0 no auxilliary pressure plot 0 

= 1 auxilliary pressure plot made. 


22 , 

22 , 

30 


XPLT.YPLT Real, length in inches of the x axis and 

C p axis respectively. 

CPI Real, maximum and minimum C n ordinates. 

CP2 p 


5. ,8. 

1 . 2 , 

- 2.0 


DTIC 


Real, Cp increments between ordinate tic marks. .4 


AMXP.ALXP, Real, experimental values of Mach number, angle .5,0., 
CLXP ,CDXP of attack, lift coefficient, drag coefficient 0.,0., 
CMXP and moment coefficient, respectively. 0. 


DCLC Real, increment in input lift coefficient, CLFT1, 0. 

required to match calculated and experimental 
lift due to difference in circulatory lift and 
integrated pressure distribution. Requires two 
runs. 
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Parameter Comments Default 

NXP,NXL Integer, total number of experimental pressure 0,0 

points and number of points on lower surface, 
respectively. If NXP=0, program will not expect 
any data. The pressure data C p (x) appears 
later on in the data set. If NXP > 0, program 
sets IPSW=1. 

KPTS,KLS Integer, total number of boundary layer data 0,0 

points and number of points on lower surface, 
respectively. If KPTS=0 program will not expect 
any data. The data <S*(x), 0(x), and Cf(x) 
follows the pressure data, if any, and appears 
later on in the data set. 

The next record in the data set is another dummy header and is read from 
the program by: 

READ( IREAD,500) 

500 FORMAT (IX) 

This record gives the names of the parameters used in controlling the numerics 
of the relaxation process, i.e., maximum number of multi-grid cycles, 
convergence criteria, ADI constants and the number of micro-meshes. It is 
followed by as many records as there are macro-meshes and they are read from 
the program as: 

READ( IREAD, 510)T0T0( I ) ,C0V0( I), P 10 (I), P 20 (I), P 30 (I), P 40(1), P 50 (I) ,GMESH( I ) 

510 F0RMAT(8F10.6) 

where: 

T0T0(I) REAL Cols 1 - 10 Maximum number of multi -grid cycles to limit 

unnecessary computations in poorly 
converging situations. For an inviscid 
calculation, 25 is a reasonable value for 
each macro-mesh. It may have to be 
increased in difficult cases to obtain the 
desired degree of convergence. For viscous 
calculations, use 2 on each macro-mesh in 
order to update boundary layers most 
frequently. 

C0V(I) REAL Cols 11-20 Convergence criterion on maximum residual .5 

x 10" 4 for each macromesh is usually 
sufficient (right justified). 
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P10(I) 

REAL 

Cols 

21 

30 

P20( I ) 

REAL 

Cols 

31 

40 

P30( I ) 

REAL 

Cols 

41 

50 

P50( I ) 

REAL 

Cols 

61 

70 

P40( I ) 

REAL 

Cols 

51 

60 


GMESH(I) REAL Cols 71-80 


These ADI parameters can be kept 
constant with a small sacrifice in 
convergence as: 


.25 

.45 

.45 


J 1.5 


Supersonic relaxation factor. For inviscid 
cases use .8 for each macro-mesh with 
incremental reduction to .5 in cases where 
divergence or oscillatory convergence is 
seen in the residual. For viscous cases, 
use .6, .7, and .8 on the three macro-meshes 
to allow for the unsettling effect of the 
boundary layer in the early stages. They 
can be reduced as convergence warrants. 


The number of micro-meshes associated with 
each macro-mesh. Since a micro-mesh is 
derived from a macro-mesh by halving the 
number of divisions in each direction, care 
must be taken to prevent the number of 
radial divisions from falling below 2. 
Similarly, the values of GMESH must be 
chosen with due regard for the values of FNC 
and FNR in that each halving must result in 
integer values for the micro-mesh divisions. 


Another dummy header identifying QC, VIS1, VIS2 and VIS3 follows read by: 


READ (IREAD,500) 
500 FORMAT (IX) 


The next record contains the values of the four parameters read by: 

READ(IREAD,510)QC,VIS1,VIS2,VIS3 
510 FORMAT (8F10.6) 


where: 

QC 


VIS1 

VIS2 


REAL Cols 1-10 Lower bound on square of local flow field 

Mach number where artificial viscosity is 
switched on. A value of .9 is suggested. 
For cases where convergence appears to hang 
up a reduction to .85 may help with little 
loss in accuracy of forces if the solution 
is second-order. 

REAL Cols 11-20 First-order artificial viscosity factor. 

Use 1.0 


REAL Cols 21-30 Switch to established order of artificial 

vi scosity 

0. = second-order 

1. = first-order 
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VIS3 REAL Cols 31-40 Second-order artificial viscosity factor. 

Use 2.0. 

Another dummy header record is read by: 

READ(IREAD,500) 

500 (IX) 

which identifies the parameters for angle of attack and Mach number variation 
that the user assigns on the next record which is read from the program as: 

READ ( IREAO , 510 ) AL1 ,AL2 .STEP ,FM1 ,FM2 ,DMACH 
510 FORMAT (8F10.6) 


AL1 

REAL 

Cols 

1-10 

Initial angle of attack in degrees. 

AL2 

REAL 

Cols 

11-20 

Final angle of attack in degrees. 

STEP 

REAL 

Cols 

21-30 

Incremental angle of attack in degrees 

FM1 

REAL 

Cols 

31-40 

Initial Mach number. 

FM2 

REAL 

Cols 

41-50 

Final Mach number. 

DMACH 

REAL 

Cols 

51-60 

Incremental Mach number (+ or -). In 


inverse mode AL1, AL2, and STEP may be left 
blank, and if only AL1 and FM1 are used, 
four other parameters are left blank. 

Similarly, another header for describing the lift coefficient variation is 
read: 


READ( IREAO, 500) 

500 FORMAT (IX) 

followed by the actual values and read from the code by: 

READ( IREAD,510)CLFT1 ,CLFT2,DCLFT 
510 FORMAT (8F10.6) 


where: 


CLFT1 

REAL 

Cols 

1-10 

Initial lift coefficient. 

CLFT2 

REAL 

Cols 

11-20 

Final lift coefficient. 

DCLFT 

REAL 

Cols 

21-30 

Incremental lift coefficient (+ or -). In 


direct mode CLFT1,CLFT2 and DCLFT may be 
left blank, and if only CLFT1 is used, two 
remaining parameters are left blank. 


Read another header by 
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READ ( IREAD , 500 ) 
500 (IX) 


Then, 

READ ( IREAD, 510 )FSYM,FNU,FNL 
510 FORMAT (8F10.6) 

where: 

FSYM REAL Cols 1-10 must be +1. or +2. 

+ symmetrical profile, only upper surface 
ordinates needed. 

- unsymmetrical profile, lower surface 
ordinates also needed. 

= +1 surface slopes are computed in program. 
= ±2 surface slopes are given. 

FNU REAL Cols 11-20 Number of points on upper surface (suggested 

minimum for definition, 30). 

FNL REAL Cols 21-30 Number of points on lower surface (suggested 

minimum for definition, 30). For 
symmetrical profile set FNL = FNU. 

Read a header which identifies items to be read by next read statement 

READ(IREAD,500) 

500 FORMAT (IX) 

READ( IREAD, 510 )TRAI L ,SL0PT 
510 FORMAT (8F10.6) 

where 

TRAIL REAL Cols 1-10 Trailing edge included angle in degrees. If 

a negative number is used, program will 
compute the actual angle. 

SLOPT REAL Cols 11-20 Trailing edge slope. If surface slopes of 

profile are given ( | FSYM | =2) , TRAIL must 
also be given and SLOPT is not required. If 
surface slopes of profile are not given, 
TRAIL will be calculated if a negative value 
is read in and SLOPT is not required; 
however, if positive TRAIL read in, SLOPT 
must also be read in. 

The next group of records defines the profile in terms of its ordinates 
and surface slopes. The surface slopes should only be prescribed if they are 
very accurately known since the program will calculate them from the 
ordinates. The x-values of the trailing edge must be equal on the upper and 
lower surface; however, the y-values need not be equal, in which case the 
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airfoil has an open trailing edge. If the profile is symmetric, only the 
cards describing the upper surface should be supplied. 

Read a header which identifies items to be read from the next FNU records 
defining the upper surface of the profile. 

READ( IREAD,500) 

500 FORMAT (IX) 

Then: 


READ(IREAD,510) X(I), Y(I), SLOP 
510 FORMAT (8F10.6) 


X(I) 

REAL 

Cols 

1-10 

X ordinates of upper surface from LE to TE. 

Y ( I ) 

REAL 

Cols 

11-20 

Y ordinates of upper surface. 

SLOP 

REAL 

Cols 

21-30 

slopes on upper surface, in degrees. Can be 
left blank unless very accurately known. 


Read a header which identifies items to be read from the next FNL records 
defining the lower surface of the profile. 

IF FSYM < 0. READ (IREAD,500) 

500 FORMAT (IX) 

Then: 


READ(IREAD,510) X(I),Y(I), SLOP 
510 FORMAT (8F10.6) 


X(I) 

REAL 

Cols 

1-10 

X ordinates of lower surface from LE to TE. 

Y ( I ) 

REAL 

Cols 

11-20 

Y ordinates of lower surface. 

SLOP 

REAL 

Cols 

21-30 

slopes on lower surface, in degrees. Can be 
left blank unless very accurately known. 


The convention for slopes is as follows: the horizontal tangent to the 

lower surface corresponds to 0°, the vertical tangent at the leading edge 
corresponds to -90°, and the horizontal tangent to the upper surface 
corresponds to -180°. 

The next two input records refer to NAMELIST statements which control the 
details of the boundary layer calculation. Default values are used for any 
unspecified parameters. These records are always included, however, for 
in viscid runs they function as dummies. 
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The first record controls seven integer parameters using the read 
statement: 


Parameter 


KT = 1 Ref. 11 
= 2 Ref. 12 
= 3 Ref. 13 


READ ( I READ , 1 1 NP ) 

Comments Default 


Transition options: 

Natural transition mode. User selects 1 

one of three possible forms. If user 

makes this choice, he must also set XTRANU, L = 0. 


= 4 Transition fixed mode. User assigns transition 

point on upper and on lower surface. 

XTRANU > 0. and XTRANL > 0. 


ICURV 
I D I V 
IH1HB 

=1 

=2 

=3 

IMCE 

=1 

=2 

ICF 

=1 

=2 


Combined mode. User sets KT=1,2, or 3 and also 
sets XTRANU > 0. and XTRANL > 0. In this mode 
program will choose the earlier occurring 
of the two transition locations, fixed or natural. 
This is a good mode of operation if transition 
locations are not known. Since the default values 
of XTRANU ,L = .10, this is default mode of code 
operation. Note: If short or long bubble leading 

edge separation is detected within the program, 
immediate transition is assumed and any user 
selected options are overridden. 


Longitudinal curvature term in turbulence 
model: = 0 term out, = 1 term in. 

Mean dilatation term in turbulence model: 

= 0 term out, = 1 term in. 

Turbulent velocity profile used with the 
lag-entrainment method (see Ref. (6)): 
relation in Green's original lag-entrainment 
method 

East, Smith and Merryman 
Green-LeBal leur 

Equilibrium entrainment model: 
relation in Green's original lag- 
entrainment method 

blend of Michel's and a free mixing layer 
Skin friction relation: 

relation in Green's original lag-entrainment 

method 

Swafford's 


1 

1 

3 


2 


2 
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Parameter Comments Default 

IOUT Output parameter for boundary layer results: -1 

=1 prints out all tables after every 

boundary layer iterate 

=0 prints out all tables after the last 

boundary layer iterate on each macro-mesh 
=-l prints out all tables after the last boundary 

layer iterate on the final macro-mesh 


The next record to be read controls 26 NAMELIST parameters, all of which 
are real, using the read statement: 



READ(IREAD,AINP) 


REYNO 

Reynolds number based on chord 

NONE 

XTRANU , 
XTRANL 

Specified transition point on upper and on 
lower surface 

.10, 

.10 

VMAX, 

VMAXT 

Sets upper bound on source velocity, V, near 
shock waves, the trailing edge, and in the wake. 

• • 

cn cn 
o o 

w 


Used to prevent divergence during early stages 
of iteration due to nonphysical boundary layer 
behavior. V = VMAX if V >VMAX. Also provides 
a reasonable model of interaction near strong shock 
waves where V could be greater than VMAX. VMAXT is 
invoked upstream and within .95 of trailing edge. 
(VSMAX is a similar parameter for the symmetrical 
source velocity downstream of the trailing edge 
and is VMAXT/ 2.) 


RQEO, Base relaxation factor applied to the .2, 
RQE1, source strength on the airfoil. One 1.0, 
RQE2 value assignable to each macro-mesh. .8 

In difficult cases drop RQE1 to .5 and 
RQE2 to .2. 

WRQO, A multiplicative factor used with RQE to 5., 
WRQ1, provide a buildup to the maximum on 1., 
WRQ2 each macro-mesh. This is particularly 1. 

important on the coarse mesh when the 
calculation starts off. In fact, for 


difficult cases, only WRQO is changed to 
4.0 or 3.5. 

RQWO, Base relaxation factor applied to the .2, 

RQW1, source strength in the wake. One value 1.0, 

RQW2 assignable to each macro-mesh. In difficult .4 

cases drop RQW1 to .4 and RQW2 to .2. 
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Parameter Comments Default 

WRQWO, A multiplicative factor used with RQW to 5., 

WRQW1, provide a buildup to the maximum on each 1., 

WRQW2 macro-mesh. Here too, the buildup is 1. 

important on the first macro-mesh and, in 
difficult cases, WRQWO is reduced to 2.5. 

RBWO, Base relaxation factor applied to the wake .2, 

RBW1 , angle. One value is assignable to each .4, 

RBW2 macro-mesh and, in difficult cases, RBW1 is .2 

dropped to .2 and RBW2 is dropped to .1. 

WRBWO, A multiplicative factor used with RBW to 5., 

WRBW1, provide a buildup to the maximum on each 1., 

WRBW2 macro-mesh. Again the buildup is important 1. 

on the first macro-mesh and, in difficult cases, 

WRBWO is reduced to 1.5. 

RTEO, Under-relaxation of a program generated 1., 

RTE1, constant used in the calculation of the 1., 

RTE2 wake source velocity is very helpful in .1, 


achieving boundary layer convergence. This 
is especially true on the finest macro-mesh and, 
in fact, no relaxation is applied to the earlier 
meshes under usual conditions. For difficult 
cases RTE2 is reduced to .05. 

At this point in the data set, experimental data may be introduced. The 
data can be of two main types. If both types are available, the pressure data 
are inserted first and are followed by the boundary layer data. For either 
type the lower surface data are read first followed by the upper surface 
data. Within these categories the data may be in any order. The user will 
have alerted the program that data is to be expected here by assigning non 
zero values to NXP and/or KPTS in the MXF NAMELIST. 

The program will then read NXP records from: 

READ(IREAD,32) XXP(I), CPXPI(I) 

32 FORMAT (2E10.7) 

where: 

XXP(I) REAL Cols 1-10 X location of pressure tap 

CPXPI(I) REAL Cols 11-20 C p at location X 

Next the program reads KPTS records from: 
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Fig. 3 Card Images of Input File 
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READ ( IREAD ,510) XDAT ( I ) ,DSDAT( I) ,THDAT ( I ) ,CF0AT( I) 
510 FORMAT (8F10.6) 


where: 


X DAT ( I ) REAL Cols 1-10 
DSDAT(I) REAL Cols 11-20 
THDAT(I) REAL Cols 21-30 
CFDAT(I) REAL Cols 31-40 


X location of data point 
5* at location X 
0 at loclation X 
Cf at location X 


If, for some reason, any of the parameters <S*, 0, or Cf are missing at some X 
location, the user should put in a dummy value of 1.0 for the parameter. The 
code will recognize this special value and not attempt to plot the point. 


When the code has processed all the information contained in the data set 
thus far, the program returns to statement number one. It is ready to process 
the data for another profile, or as many as computing time will allow. At 
this point the user turns the program off by adding the first three records of 
a final data file. On the first record, which specifies the new title, the 
user (at his option) could type the phrase, "End of Calculation". This is 
appropriate since this phrase will be the last thing printed by the code. The 
next record, a dummy labeled FNC is followed by the third and final record 
with the value 0. for FNC. The program uses this value of FNC to branch to a 
STOP. 


3.3 Output Description 

Much of the output for either an inviscid or viscous calculation is 
controlled by the output parameter IOPUT as indicated in Table I. A 
reasonable amount of paper is generated when this parameter equals 1. For the 
boundary layer characteristics the output parameter is IOUT, and if -1 is 
used, this material will only be printed for the last boundary layer 
calculation on the last macro-mesh. The included sample case was run with 
IOPUT = 1 and IOUT = -1. The first item in the Appendix is the data file 
which produced the output given as the second item in the Appendix. 

The following description of the output makes reference to the sample 
output included in the Appendix. The solution output is provided in a series 
of tables and graphs. Since the graphical output was described in some detail 
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in Section 2.2, we will only describe the tabular output. The tabular results 
are presented in the following eight major blocks: 

1. Program identification block, 

2. Title or profile description, 

3. Initial data, 

4. Coarse grid convergence data, 

5. Boundary layer solution results, 

6. Inviscid iteration quantities, 

7. Inviscid solution quantities 

8. Fine grid convergence data 

The position of each of these blocks is indicated by the corresponding numbers 
in the sample output. 

3.3.1 Program Identification (1) 

This block identifies the version of the code being used and the date 
that it was released. 

3.3.2 Title or Profile Description (2) 

The second block of data, (2), identifies the profile and case being run, 
the prescribed lift coefficient (or angle of attack if a direct problem were 
being run) and the free stream Mach number. 

3.3.3 Initial Data (3) 

This set of data contains all the additional information required to run 
the program. The data in Block (3) includes tables giving the airfoil 
geometry, values of parameters controlling the inviscid computation, and the 
viscous calculations. It also includes a tabular listing of experimental 
data, if provided by the user, which is used for comparison with the numerical 
solution in the graphical output. The airfoil geometry table in this block 
provides a listing of arc length, normalized (by chord length) coordinates, 
surface slope, and curvature. The table is followed by the values of four 
important geometric properties of the airfoil, namely: leading edge radius, 

trailing edge angle, leading edge slope and trailing edge slope. 



Diagnostic checks are made at selected points in the program to insure 
its proper operation. If preset conditions are not met, the program will stop 
and a diagnostic message will be printed. The first such message will appear 
following the airfoil identifier (block 2) if, for a viscous calculation, the 
trailing edge values of X are unequal. Another program stop will occur and a 
diagnostic message will be written following the airfoil geometry table, if 
more than 50 iterations are needed to converge the circle plane mapping. In 
this instance the user is urged to check the profile ordinates for a 
mistype. Scanning the profile curvature in the table can be helpful in 
locating a bad point. 

The inviscid calculation parameters appear in the mext set of entries 
following the geometry table. These include: mode of code operation (LMODE, 

NINV.and IRH ) ; the artificial viscosity parameters (QC, VI SI , VIS2, and VIS3); 
the output parameter, IOPUT, which controls the inviscid solution output; the 
multi -grid-ADI parameters. 

The parameters controlling the viscous calculation are given in the next 
set of entries. Subheadings in this data block are self-explanatory. 

The experimental data for pressure coefficient that is supplied for 
plotting purposes is listed in a table of X and Cp. Since AMXP and FM1 are 
equal (see data file), there is no Mach number correction to the data, and the 
third column is the same as the second, and the corrected and uncorrected 
forces are also equal. 

The scaling parameters that control an auxiliary plot showing a 
comparison of experimental and theoretical surface pressures are listed 
following the pressure coefficient table. The fact that experimental Cp data 
is given triggers the auxiliary plot. The scalings given in the sample output 
are default values. 

If experimental data for boundary layer properties are supplied (for 
incorporation into the graphical output) a table of such data will be included 
in the output following the pressure data. Since such data was supplied in 
the sample test case, a table of X, 6*, 6, and Cf appears at the end of the 
block (3) ouptut. 
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3.3.4 Coarse Grid Convergence Data (4) 

The program is set-up to use three macro-grids in carrying out the 
solution. The convergence history of the solution on the first two grids is 
given in the fourth (4) block of output data. The convergence tables 
summarize how various flow quantities change as the iterative solution 
progresses . 

In an uneventful situation the only output for the first macro-mesh is 
the convergence summary table (output section 4). However, before discussing 
this table, we offer a few words about less routine operation. At various 
points in the program several parameters are checked to see if they fall 
within reasonable physical limits. If these limits are exceeded, a limit 
value is assigned and a flag is printed. These flags may appear on any mesh 
early-on and just be transient in that they disappear as convergence 
proceeds. If they persist, however, the fix is to reduce the boundary layer 
relaxation factors on the problem macro-mesh. 

A more frequent occurrence is, as shown in the sample case, the 
indication of a movement of the transition point. This can be due to a change 
in the solution or due to the users specified transition point being too close 
to the stagnation point. If a shift occurs, the code smooths the source 
strength Q in the vicinity of the transition point to soften its effect on the 
solution. The relevant information as well as the unsmoothed and smoothed 
values of Q and their locations are also printed. Under some conditions the 
transition point moves back and forth and the solution doesn't converge. 
Usually this problem can be brought under control by reducing RQE and/or using 
KT = 4. 

In the convergence tables the quantity NIT is the number of complete 
multigrid cycles; RMAX is the maximum residual; VTOP, VBOT, and VSYM are 
trailing edge values of the source velocities on the top and bottom surfaces 
and the symmetrical wake value; DELT and DELB are the boundary layer thickness 
for the top and bottom surfaces at the trailing edge or at SEPT and SEPB 
points upstream of the trailing edge in separated flow; DUT and DUB show 
respectively the maximum difference between UEI and UEV on the top and bottom 
side of the airfoil and wake. The integers listed under I-J are the grid 
labels identifying the grid point on the airfoil or wake where the maximum 
difference DUT or DUB occurs. An unsigned integer indicates that the maximum 
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occurs on the airfoil at the point I noted in the table. A minus sign 
indicates it occurs in the wake at the point J in the table. The values of I 
fall in the range 3 < I < NC + 3 where 1=3 and NC + 3 are the lower and 
upper surface trailing edge points respectively. The index J falls in the 
range 2 < J < NR + 2 where J = 2 is the wake point at infinity and J = NR + 2 
is the wake point at the trailing edge. The quantity, BOMX is the maximum 
difference between boundary layer iterates of the wake angle which is located 
at grid point, JMX, on the wake axis. Note that the code generates an initial 
approximation to the source strength Q to obtain the first boundary layer 
solution. It then updates this solution fifteen successive times before 
returning the updated boundary conditions to the inviscid calculation, and 
thus, ITBL starts at 15. Note also that the last line (22) in the table, 
ITBL=36, which corresponds the maximum number of iterates specified for this 
macro-mesh, was the cause for ending the first mesh calculation since the 
convergence checks on RMAX, AL, CPTE, DUT, and DUB have not been met. 

However, since the convergence is generally monatonic, this is an acceptable 
first mesh solution. 

Mesh division is undergone, interpolation of pertinent arrays 
accomplished, a solution obtained, and the second macro-mesh results are 
printed.. At mesh division it is possible for the transition point to shift 
simply because there are more nodal points on the refined mesh as shown in the 
sample case. The convergence table shows that almost all the convergence 
conditions have been met in 26 boundary layer iterations; After the first 
mesh division, five successive boundary layer passes are made before returning 
updated boundary conditions to the inviscid calculation, and, thus, ITBL 
starts out at 5. This is done to help the solution negotiate mesh division. 

On the second mesh division only two successive boundary layer passes are 
made. 

3.3.5 Boundary Layer Solution Results (5). 

On the third mesh the full boundary layer output is presented. The first 
table lists the upper surface boundary layer parameters. H is Green's shape 
parameter and CE is the entrainment coefficient. Printed at the bottom of the 
table are the location and value of the maximum difference in UEI and UEV as 
well as the average difference. 
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A LAMBDA flag is printed when the effects of the ICURV and IDIV 
turbulence modeling terms exceed reasonable physical values. Next is shown 
the iterative solutions for tt and 6 at the trailing edge given the values of 
6* and Cf from the boundary layer solution. A parallel set of results for the 
lower surface follows. Next, a short table of various drag quantities is 
presented. The quantity, CDINF is the profile drag coefficient calculated 
from momentum thickness far downstream; CDSY is the profile drag coefficient 
calculated from Squire-Young formula based on trailing edge momentum 
thickness; CDINFT and CDIFB are the top and bottom surface contributions to 
CDINF; CDF is the total skin friction drag coefficient while CDFT and CDFB are 
the top and bottom surface contributions to CDF. 

3.3.6 Inviscid Iteration Quantities (6) 

The next table, identified as iterative solution, lists the ADI and 
artificial viscosity parameters, and the mesh geometry for the finest macro- 
mesh. This is followed by the relaxation solution for the last boundary layer 
iterate. The table headings are mainly self-explanatory. The tenth column is 
the circulation or angle of attack for LM0DE=0 or -1, respectively. The next 
column, headed work, is the cumulative value of the equivalent SLOR steps that 
have taken place from the first relaxation step on the current macro-mesh to 
the current cycle. The last column refers to the number of supersonic points 
in the flow field. The next two lines specify: the average residual from the 

first relaxation step on the current macro-mesh; the average residual from the 
last relaxation step on the current macro-mesh; the cumulative work, an 
overall convergence rate based on the cumulative work; and an overall 
convergence rate based on the cumulative number of cycles. Another line is 
only printed when the Hafez entropy correction is made, IRH=1. It identifies 
the correction at the body surface for the last shock point on the surface. 

The index I refers to the location of the point where the correction is made; 
EDSR is equal to the exponential of the entropy jump across the shock wave 
divided by the gas constant; the next two terms are the Mach numbers upstream 
and downstream of the shock respectively. If a shock were present on the 
lower surface, this line would be printed out first and a second line would be 
printed if a shock were on the upper surface. 

The next line gives the macro-mesh number, the cumulative number of 
boundary layer iterates and the cumulative number of multi grid cycles for this 
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mesh. Following this is a series of lines giving the approximate shock 
location, and pressure computed from the "exact" and irrotational jump 
conditions based on the Mach number at the first supersonic point upstream of 
the shock wave at the body surface. These four lines are printed for every 
shock on the profile. 

3.3.7 Inviscid Solution Quantities (7) 

The first two lines in this block provide: the angle of attack 

(calculated in the sample case), Mach No., CL calculated from the circulation 
(given in the sample case). The next line gives the CL, CD, CDB (base drag 
included), and moment -computed for the composite surface pressure 
distribution. 

The next table summarizes the solution for the airfoil surface properties 
in terms of the coordinates of the nodal points on the final mesh starting at 
the lower surface trailing edge and moving around the profile to the upper 
surface trailing edge. In addition to the arc length, it lists the 
uncorrected or inviscid: Mach No., velocity ratio, density, mass flow, and 

pressure coefficient. The last column is the corrected or composite pressure 
distribution. 

3.3.8 Fine Grid Convergence Data (8) 

The final mesh convergence summary table is given in this block. The 
nomenclature is the same as used in the coarse grid convergence tables 
discussed in subsection 3.3.4. The Richardson extrapolation of the angle of 
attack, CL, CDB and CM from the last two macro-meshes is the final calculated 
result to be printed. 

The formal end of the calculation repeats the program title block and 
prints the end of calculation title. 

Both the printed and graphical output were generated with the Xerox 9700 
printer. In the case of the graphics, rotation of the plots and their scaling 
reflect changes made by the 9700 and not contained within the program. 
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APPENDIX 
SAMPLE CASE 


DATA FILE 


RAE 2822 CASE 9 SAMPLE CASE FOR LANGLEY FINAL REPORT 

FNC FNR FMESH 

40. 8. 3. 

11-1 1 
&MXF DCLC= . 0053 , 

AMXP- .730,ALXP=3 . 1 9 , CLXP= . 803 , CDXP= . 01 68 , CMXP=- . 099 , 

NXP= 102 , NXL = 50 , KPTS= 1 1 , KLS = 2 , &END 

TOTO COVO P10 P20 P30 P40 P50 GMESH 


2. 

. 5E-4 

.25 

.45 

.45 

.6 

2. 

. 5E-4 

.25 

.45 

.45 

.7 

2. 

. 5E-4 

.25 

.45 

.45 

.8 

QC 

VIS1 

VIS2 

VI S3 



.90 

1 . 

0.00 

2.00 



AL 1 

AL2 

STEP 

FM1 

FM2 

DMACH 

2.23 

.000 

.000 

.730 

.000 

.000 

CLFT1 

CLFT2 

DCLFT 




.803 

.000 

.000 




FSYM 

FNU 

FNL 




-1 . 

65. 

65. 




TRAIL 

SLOPT 






-1.5 

X Y UPPER SURFACE 


0.000000 

0 . 000000 

0.000602 

0.003165 

0.002408 

0.006306 

0.005412 

0.009416 

0.009607 

0.012480 

0.014984 

0.015489 

0.021530 

0.018441 

0.029228 

0.021348 

0.038060 

0.024219 

0.048005 

0.027062 

0.059039 

0.029874 

0.071136 

0.032644 

0.084265 

0.035360 

0.098396 

0.03801 1 

0. 113495 

0.040585 

0. 129524 

0.043071 

0. 146447 

0.045457 

0. 164221 

0.047729 

0. 182803 

0.049874 

0.202150 

0.051885 

0.222215 

0.053753 

0.242949 

0.055470 

0.264302 

0.057026 

0.286222 

0.C58414 

0.308658 

0.059629 

0.331555 

0.060660 

0.354858 

0.061497 

0.378510 

0.062133 

0.402455 

0.062562 

0.426635 

0.062779 

0.450991 

0.062774 

0.475466 

0.062530 

0 . 500000 

0.062029 

0.524534 

0.061254 

0.549009 

0.060194 

0.573365 

0.058845 

0.597545 

0.057218 

0.621490 

0.055344 

0.645142 

0.053258 

0.668445 

0.050993 

0.691342 

0.048575 

0.713778 

0.046029 

0.735698 

0.043377 

0.757051 

0.040641 

0.777785 

0.037847 

0.797850 

0.035017 

0.817197 

0.032176 

0.835779 

0.029347 

0.853553 

0.026554 

0.870476 

0.023817 

0.886505 

0.021153 


48 



0.001604 

0.018580 

0.915735 

0.0161 13 

0.928864 

0.013769 

0.940961 

0.01 1562 

0.951995 

0.009508 

0.961940 

0.007622 

0.970772 

0.005915 

0.978470 

0.004401 

0.985016 

0.003092 

0.990393 

0.002001 

0.994588 

0.001 137 

0.997592 

0.000510 

0.999398 

0.000128 

1 .000000 

0 . 000000 

X 

Y 

0.000000 

0.000000 

0 . 000602 

- 0.003160 

0.002408 

- 0.006308 

0.005412 

- 0.009443 

0.009607 

- 0.012559 

0.014984 

- 0.015649 

0.021530 

- 0.018707 

0.029228 

- 0.021722 

0.038060 

- 0.024685 

0.048005 

- 0.027586 

0.059039 

- 0.030416 

0.071136 

- 0.033170 

0.084265 

- 0.035843 

0.098396 

- 0.038431 

0 . 113495 

- 0.040929 

0 . 129524 

- 0.043326 

0 . 146447 

- 0.045610 

O , 164221 

- 0.047773 

0. 182803 

- 0.049805 

0.202150 

- 0.051694 

0.222215 

- 0.053427 

0.242949 

- 0.054994 

0.264302 

- 0.056376 

0.286222 

- 0.057547 

0.308658 

- 0.058459 

0.331555 

- 0.059046 

0.354858 

- 0.059236 

0.378510 

- 0.058974 

0.402455 

- 0.058224 

0.426635 

- 0.056979 

0.450991 

- 0.055257 

0.475466 

- 0.053099 

0 . 500000 

- 0.050563 

0.524534 

- 0.047719 

0.549009 

- 0.044642 

0.573365 

- 0.041397 

0.597545 

- 0.038043 

0.621490 

- 0.034631 

0.645142 

- 0.031207 

0.668445 

- 0.027814 

0.691342 

- 0.024495 

0.713778 

- 0.021289 

0.735698 

- 0.018232 

0.757051 

- 0.015357 

0.777785 

- 0.012690 

0.797850 

- 0.010244 

0.817197 

- 0.008027 

0.835779 

- 0.006048 

0.853553 

- 0.004314 

0.870476 

- 0.002829 

0.886505 

- 0.001592 

0.901604 

- 0.000600 

0.915735 

0.000157 

0.928864 

0.000694 

0.940961 

0.001033 

0.951995 

0.001 197 

0.961940 

0.001212 

0.970772 

0.001 1 12 

0.978470 

0.000935 

0.985016 

0.000719 


LOWER 


SURFACE 
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0.990393 0.000497 

0.994588 0.000296 

0.997592 0.000137 

0.999398 0.000035 

1 . 000000 0 . 000000 
&IINP SEND 

&AINP REYN0=6.5E6,XTRANU 
WRQ0=4 . , WRQW0=3 . , WRBW0=2 


030 , XTRANL= . 030, 


ROE 1= . 60, RQW1= .40. RBW1= . 20, 
RQE2= . 3 , RQW2= . 20, RBW2= . 10, SEND 


.0000 

0.7963 

.0002 

1 .0157 

.0007 

1.0912 

.0015 

1.1315 

.0025 

1 . 1257 

.0036 

1 ,0944 

.0047 

1 .0467 

.0060 

0.9971 

.0073 

0.9400 

.0087 

0.8967 

.0104 

0.8441 

.0125 

0.7849 

.0146 

0.7309 

.0186 

0.6501 

.0271 

0.5652 

.0375 

0.4393 

.0500 

0.3447 

.0625 

0.2776 

.0750 

0.2263 

. 1500 

0.0203 

.2000 

-0.0641 

.2500 

-0. 1462 

.3000 

-0.2246 

.3250 

-0.2635 

. 3500 

-0.2960 

. 3750 

-0.3162 

.4000 

-0.3187 

.4250 

-0.2987 

.4 500 

-0.2673 

,4750 

-0.2292 

,5000 

-0. 1836 

, 5250 

-0. 1441 

,5500 

-0.0971 

,5750 

-0.0522 

,6000 

-0.0108 

,6196 

0.0168 

,6500 

0.0594 

,7000 

0. 1200 

7500 

0. 1786 

7750 

0.2038 

8000 

0.2281 

8250 

0.2471 

8500 

0.2657 

8750 

0.2792 

9000 

0.2917 

9250 

0.2965 

9500 

0.2922 

9750 

0.2640 

9875 

0.2381 

9937 

0.2100 

0002 

0.3761 

0007 

0.0491 

0015 

-0. 1075 

0025 

-0.3348 

0036 

-0.4609 

0047 

-0.5852 

0060 

-0.6827 

0073 

-0.7675 

0087 

-0.8198 

0104 

-0.8572 

0125 

-C.91 13 

0146 

-0.9608 

0187 

- 1 .0439 

0271 

-0.9157 

037 5 

-1.1741 
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.0500 -1.1973 
.0625 -1.1771 
.0750 -1.1697 
. 1000 -1.1 123 
.1500 -1.0981 
.2000 -1.0834 
.2208 -1.0890 
.2500 -1.1110 
.3000 -1.1535 
.3250 -1 . 1658 
.3500 -1.1727 
.3750 -1.1917 
.4000 -1.2085 
.4250 -1.2290 
.4500 -1.2420 
.4750 -1.2593 
.5000 -1.2853 
.5250 -1.3097 
.5500 -1.1249 
.5750 -0.7552 
.6000 -0.6338 
.6196 -0,5681 
.6500 -0.4949 
.6771 -0.4370 
.7000 -0.3915 
.7500 -0.2983 
.7750 -0.2535 
.8000 -0.2103 
.8250 -0.1654 
.8500 -0.1226 
.8750 -0.0760 
.9000 -0.0325 


.9250 

0.0114 



.9500 

0.0570 



.9750 

0.1012 



.9875 

0. 1232 



.9937 

0. 1347 



. 152 

.000325 

.000196 

,003983 

1 .00 

.003233 

.001982 

0 . 

. 179 

.000303 

.000378 

.002700 

.319 

.001060 

.000516 

.002723 

.404 

.001311 

.000626 

.002489 

,498 

,001533 

.000731 

.002533 

.574 

.003395 

.001418 

.001049 

.65 

.004956 

.002043 

.000901 

. 75 

.006884 

.002999 

.000831 

.90 

.008958 

.004360 

.000959 

1 .00 

.014433 

.006510 

0 . 

END OF 

CALCULATION 




FNC 

0 . 
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SAMPLE OUTPUT 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxmxxxxxxxxxxxmx 

X X 

X GRUMMAK AEROSPACE CORP. PROGRAM GRUMFOIL (MCMJ-9N ) 2-22-85 X 
X X 

X AIRFOIL ANALYSIS IN TRANSONIC VISCOUS FLOW X 

X X 

X DEVELOPED BY R. MELNIK, R.CHOH, H. MEAD, AND A. JAMESON X 

X X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

RAE 2822 CASE 9 SAMPLE CASE FOR LANGLEY FINAL REPORT 
INITIAL DATA CLC- 0.8030 MACH NO-. 730 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

AIRFOIL COORDINATES 


ARC LENGTH 

X 

Y 

THETA 

KURV 

-1.0157 

1.0000 

0.0 

-3.3511 

-1.5337 

-1.0151 

0.9994 

0.0000 

-3.3042 

-1.5337 

-1.0133 

0.9976 

0.0001 

-3.1570 

-1.5337 

-1.0103 

0.9946 

0.0003 

-2.9077 

-1.3588 

-1.0061 

0.9904 

0.0005 

-2.5777 

-1.3843 

-1.0007 

0.9850 

0.0007 

-2.1508 

-1.3846 

-0.9941 

0.9785 

0.0009 

-1.6280 

-1.4018 

-0.9B64 

0.9708 

0.0011 

-1.0048 

-1.4235 

-0.9776 

0.9619 

0.0012 

-0.2968 

-1.3746 

-0.9677 

0.9520 

0.0012 

0.4611 

-1.2855 

-0.9566 

0.9410 

0.0010 

1.2259 

-1.1337 

-0.9445 

0,9289 

0.0007 

1.9707 

-1.0145 

-0.9314 

0.9157 

0.0002 

2.7033 

-0.9318 

-0.9172 

0.9016 

-0.0006 

3.4148 

-0.8231 

-0.9021 

0.8865 

-0.0016 

4.0903 

-0.7353 

-0.8860 

0.8705 

-0.0028 

4.7195 

-0.6307 

-0.8690 

0.8536 

-0.0043 

5.2988 

-0.5597 

-0.8512 

0. 8358 

-0.0060 

5.8322 

-0.4829 

-0.8325 

0.8172 

-0.0080 

6.3145 

-0.4182 

-0.8130 

0.7979 

-0.0102 

6.7489 

-0.3604 

-0.7928 

0.7778 

-0.0127 

7.1444 

-0.3226 

-0.7719 

0.7571 

-0.0154 

7.5071 

-0.2831 

-0.7504 

0.7357 

-0.0182 

7.8159 

-0.2173 

-0.7282 

0.7138 

-0.0213 

8.0465 

-0.1495 

-0.7056 

0.6913 

-0.0245 

8.2027 

-0.0880 

-0.6824 

0.6684 

-0.0278 

8.2796 

-0.0281 

-0.6589 

0.6451 

-0.0312 

8.2745 

0.0357 

-0.6350 

0.6215 

-0.0346 

8.1872 

0.0919 

-0.6108 

0.5975 

-0.0380 

8.0183 

0.1516 

-0.5864 

0.5734 

-0.0414 

7.7608 

0.2165 

-0.561B 

0.5490 

-0.0446 

7.3976 

0.2994 

-0.5371 

0.5245 

-0.0477 

6.9132 

0.3861 

-0.5124 

0.5000 

-0.0506 

6.2835 

0.5040 

-0.4878 

0.4755 

-0.0531 

5.4938 

0.6136 

-0.4632 

0.4510 

-0.0553 

4.5609 

0.7117 

-0.4388 

0.4266 

-0.0570 

3.5085 

0.7928 

-0.4146 

0.4025 

-0.0582 

2.3753 

0.8410 

-0.3906 

0.3785 

-0.0590 

1.2086 

0.8591 

-0.3670 

0.3549 

-0.0592 

0.0690 

0.8226 

-0.3437 

0.3316 

-0.0590 

-0.9903 

0.7642 

-0.3208 

0.3087 

-0.05B5 

-1.9236 

0.6581 

-0.2983 

0.2862 

-0.0575 

-2.712B 

0.5688 

-0.2764 

0.2643 

-0.0564 

-3.3905 

0.5090 

-0.2550 

0.2429 

-0.0550 

-4.0162 

0.5117 

-0.2342 

0.2222 

-0.0534 

-4.6287 

0.5164 

-0.2140 

0.2021 

-0.0517 

-5.2538 

0.5670 

-0.1946 

0.1828 

-0.0498 

-5.9064 

0.6049 

-0.1759 

0.1642 

-0.0478 

-6.5853 

0.6628 

-0.1580 

0.1464 

-0.0456 

-7.3049 

0.7401 

-0.1409 

0.1295 

-0.0433 

-8.0875 

0.8598 

-0.1247 

0.1135 

-0.0409 

-8.9409 

0.9782 

-0.1094 

0.0984 

-0.0384 

-9.8713 

1.1438 

-0.0950 

0.0843 

-0.0358 

-10.9217 

1.4085 

-0.0816 

0.0711 

-0.0332 

-12.1396 

1.7646 

-0.0692 

0.0590 

-0.0304 

-13.5697 

2.2593 

-0.0578 

0.0480 

-0.0276 

-15.2787 

2.9778 

-0.0475 

0.0381 

-0.0247 

-17.3455 

3.9866 

-0.0382 

0.0292 

-0.0217 

-19.8826 

5.5206 

-0.0299 

0.0215 

-0.0187 

-23.0924 

8.0326 

-0.0227 

0.0150 

-0.0156 

-27.2729 

12.1662 

-0.0165 

0.0096 

-0.0126 

-32.9483 

19.7790 

-0.0112 

0.0054 

-0.0094 

-40.9610 

33.7494 

-0.0069 

0.0024 

-0.0063 

-52.5599 

59.4550 

-0.0032 

0.0006 

-0.0032 

-69.2299 

101.2577 

0.0 

0.0 

0.0 

-89.9793 

126.0630 

0.0032 

0.0006 

0.0032 

-110.7857 

101.5355 

0.0069 

0.0024 

0.0063 

-127.5640 

60.5029 

0.0112 

0.0054 

0.0094 

-139.3802 

34.8578 

0.0164 

0.0096 

0.0125 

-147.6151 

20.4829 

0.0226 

0.0150 

0.0155 

-153.4577 

12.6174 

0.0298 

0.0215 

0.0184 

-157.6677 

7.8488 

0.0380 

0.0292 

0.0213 

-160.7500 

5.2272 

0.0473 

0.0381 

0.0242 

-163.0861 

3.5538 

0.0576 

0.0480 

0.0271 

-164.9134 

2.6128 

0.0690 

0.0590 

0.0299 

-166.4288 

2.0329 

0.0814 

0.0711 

0.0326 

-167.7276 

1.6207 

0.0948 

0.0843 

0.0354 

-168.8585 

1.3236 

0.1092 

0.0984 

0.0380 

-169.8620 

1.1129 

0.1245 

0.1135 

0.0406 

-170.7642 

0.9430 

0.1407 

0.1295 

0.0431 

-171.5855 

0.8246 

0.1578 

0.1464 

0.0455 

-172.3493 

0.7354 

0.1757 

0.1642 

0.0477 

-173.0706 

0.6697 

0.1944 

0.1828 

0.0499 

-173.7457 

0.5902 

0.2139 

0.2021 

0.0519 

-174.3771 

0.5429 

0.2341 

0.2222 

0.0538 

-174.9762 

0.4950 

0.2549 

0.2429 

0.0555 

-175.5508 

0.4689 

0.2763 

0.2643 

0.0570 

-176.1075 

0.4387 

0.2982 

0.2862 

0.0584 

-176.6395 

0.4069 

0.3207 

0.3087 

0.0596 

-177.1598 

0.4016 

0.3436 

0.3316 

0.0607 

-177.6822 

0.3939 

0.3669 

0.3549 

0.0615 

-178.2013 

0.3832 

0.3906 

0.3785 

0.0621 

-178.7167 

0.3771 

0.4145 

0.4025 

0.0626 

-179.2287 

0.3692 

0.4387 

0.4266 

0.0628 

-179.7443 

0.3751 

0.4631- 

0.4510 

0.0628 

-180.2849 

0.3996 

0.4876 

0.4755 

0.0625 

-180.8635 

0.4255 


52 


0.5121 

0.5000 

0.0620 

-181.4634 

0.4563 

0.5366 

0.5245 

0.0613 

-182.1403 

0.4778 

0.5611 

0.5490 

0.0602 

-182.8239 

0.4962 

0.5855 

0.5734 

0.0588 

-183.5161 

0.4943 

0.6098 

0.5975 

0.0572 

-184.1730 

0.4519 

0.6338 

0.6215 

0.0553 

-184.7680 

0.4128 

C.6575 

0.6451 

C . 0533 

-185.3C37 

0 . 3748 

0.6809 

0.6684 

0.0510 

-185.7959 

0.3589 

0.7040 

C . 6913 

0.0486 

-186.2564 

0.3393 

0.7266 

0.7138 

0.0460 

-186.6901 

0.3312 

0.7486 

0.7357 

0.0434 

-187.1055 

0.3255 

0.7702 

0.7571 

0.0406 

-187.4931 

0.3029 

0.7911 

0.7778 

0.0378 

-1B7.8562 

0.3029 

0.8113 

0.7979 

0.0350 

-188.1963 

0.2830 

0.8309 

0.8172 

0.0322 

-188.5105 

0.2780 

0.8497 

0.8358 

0.0293 

-188.7987 

0.2572 

0.8677 

0.8536 

0.0266 

-189.0612 

0.2521 

0.8848 

0.8705 

0.0238 

-189.3149 

0.2644 

0.9011 

0.8865 

0.0212 

-189.5559 

0.2534 

0.9164 

0.9016 

0.0186 

-189.7894 

0.2787 

0.9307 

0.9157 

0.0161 

-190.0155 

0.2717 

0.9441 

0.9289 

0.0138 

-190.2333 

0.2983 

0.9564 

0.9410 

0.0116 

-190.4467 

0.3075 

0.9676 

0.9520 

0.0095 

-190.6429 

0.3027 

0.9777 

0.9619 

0.0076 

-190.8412 

0.3812 

0.9867 

0.9708 

0.0059 

-191.0362 

0.3754 

0.9946 

0.9785 

0.0044 

-191.2231 

0.4561 

1.0012 

0.9850 

0.0031 

-191.3914 

0.4239 

1.0067 

0.9904 

0.0020 

-191.5597 

0.6469 

1.0110 

0.9946 

0.0011 

-191.7149 

0.6181 

1.0141 

0.9976 

0.0005 

-191.8847 

1.3135 

1.0159 

0.9994 

0.0001 

-191.9901 

1.3135 

1.0165 

1.0000 

0.0 

-192.0189 

1.3135 

LE RADIUS 

TE ANGLE 

LE SLOPE 

TE SLOPE 


0.0080 

8.6678 

0.0004 

-0.1349 



INVISCID CALCULATION PARAMETERS 

LHODE— 1 NINV- 4 IRH- 0 QC-0.90 VIS1-1.00 VIS2-0.C 


MULTIPLE 
GRIDS COV 

3 .50D-04 

4 .50D-04 

5 .500-04 


ADI FACTORS 


0.25 0.45 


P3 P4 
0.45 0.60 

0.45 0.70 

0.45 0.80 


VISCID CALCULATION PARAMETERS RE-.6500D+07 IOUT— 1 

MAX NO. BL ITERATIONS PER MESH ITBLH0- 22 ITBLM1* 22 ITBLM2- 30 

TRANSITION PARAMETERS 
XTRANU-0.030 XTRANL-0.C30 KT-1 

B.L. MODELLING PARAMETERS 

IMCE-2 IH1HB-3 ICURV-1 IDIV-1 ICF-2 

CONVERGENCE CRITERIA 

COVA-.005 COVL-.005 COVP-.001 DUCR«.010 

RELAXATION PARAMETERS 

HROg - ROE0-0.200 HR01 -1.00 ROE1-0.600 HRQ2 -1.00 ROE2-0.300 

HROH0- 3.00 RQW0-0.200 HROH1-1.00 RQH1-0.400 HRQW2-1.00 RQW2-0.200 

HRBN0- 2.00 RBH0-0.200 HRBH1-1.00 RBH1-0.200 HRBH2-1.00 RBH2-0.100 

RTE0-1.000 RTE1-1.000 RTE2-0.100 

UPPER BOUNDS 

VMAX-0.50 VMAXT-0.50 VSMAX-0.25 


PLOTTING EXPERIMENTAL PRESSURES WITH 
MACH NO. CORRECTION 


CP-EXP-COR 

0.7963 

1.0157 

1.0912 

1.1315 

1.1257 

1.0944 

1.0467 

0.9971 

0.9400 

0.8967 

0.8441 

0.7849 

0.7309 

0.5501 

0.5652 

0.4393 

0.3447 

0.2776 

0.2263 

0.0203 
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xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


MESH NUMBER 1 RESULTS 


LOWER SURFACE TRANSITION POINT HAS HOVED FROM I- 19 TO I- 20 
SOURCE STRENGTH HAS BEEN SMOOTHED ISTAG* 22.5 ILE- 23 


ITBL- 16 


I X 0 

20 0.4335D-01 0.4970D-04 
19 0.7968D-01 0.6677D-04 
18 0. I256D+00 0.285OD-03 
17 0. 1799D+00 0.3977D-03 


I X OSMTH 

20 0.4335D-01 0.4970D-04 
19 0.7968D-01 0.1674D-03 
18 0.1256D+00 C.2850D-03 
17 0. 1799D+00 0.3977D-03 


CONVERGENCE SUMMARY TABLE MESH 1 


ITBL 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 


NIT RMAX 
4 0.77D-02 
6 0.10D-01 

8 0.9BD-02 

10 0.59D-02 
12 0.16D-02 

14 -0. 14D-02 
16 -0.18D-02 
18 -0.22D-02 
20 -C.24D-G2 
22 -0.21D-02 
24 “0. 14D-02 
26 - 0 . 100-02 
28 -0.90D-03 
30 -0.87D-03 
32 — 0 . 8CD-03 
34 -0.73D-03 
" 0.67D-03 

•O.60D-03 
■0.53D-03 
•0.68D-03 
0.64D-03 
0.19D-03 


36 

38 

40 

42 

44 

46 



DELB SEPT SEPB 


0 


DUT I-J 
.031 41 

.186 32 

.307 33 
.325 34 
.282 34 

.208 34 

.122 34 

.083 36 

.092 36 

.087 36 

.075 36 
.065 36 

.059 36 

.054 36 
.050 36 

.045 36 
.041 36 

.037 36 

.033 36 

.029 36 
.027 36 


DUB I-J 
.016 -8 
.279 -4 

.267 -4 

.124 13 

.038 18 

.081 -3 

.088 -3 
.040 16 

.040 -4 

.035 -4 

.023 -4 
.016 18 
.014 18 

.013 18 

.011 18 
.010 18 
.009 18 

.008 18 
.008 18 
.007 -4 

.006 -4 


0 .023 36 .007 


BJMX JHX 
0.0000 ~ 
-.0071 
-.0078 
-.0067 
-.0054 
-.0040 
0.0028 
0.0031 
0.0036 
0.0039 
0.0036 
0.0030 
0.0027 
0.0024 
0.0022 
0.0021 
0.0019 
0.0018 
0.0017 
0.0017 
0.0016 
0.0015 


MESH NUMBER 2 RESULTS 


LOWER SURFACE TRANSITION POINT HAS MOVED FROM I« 37 TO I* 36 
SOURCE STRENGTH HAS BEEN SMOOTHED ISTAG- 41.5 ILE- 42 


ITBL- 


I X Q 

36 0.6018D-01 0. 5843D-04 
35 0. 7962D-01 0.2029D-C3 
34 0. 1015D+00 0.2296D-03 
33 0. 1256D+00 0.2563D-03 


I X OSMTH 

36 0.601BD-01 0.5843D-04 
35 0.7962D-01 0.1431D-03 
34 0. 1015D+00 0.2089D-03 
33 0. 1256D+00 0. 2563D-03 


USER’S LOWER SURFACE TRANSITION POINT TOO CLOSE TO STAGNATION POINT 
RESET TO XTL-0.060 FOR THIS MESH 


CONVERGENCE SUMMARY TABLE MESH 2 


ITBL 

NIT 

RMAX 

AL/CL 

CDB 

CPTE 

VTOP 

VBOT VSYM 

DELT DELB SEPT 

SEPB DUT 

T-J 

DUB 

I-J 

BJMX JHX 

5 

4 

-0.74D-02 

3.069 

0.0164 

0.3741 

0.459D-01 

0.211D-01 -0.253D+00 

.03343 .01329 

0 

0 .142 

67 

.049 

-16 

0.0041 

12 

6 

6 

-0.36D-02 

2.921 

0.0175 

0.2454 

0.111D-01 

0. 154D-01 -0. 106D+00 

.03376 .01311 

2 

0 .177 

66 

.059 

-14 

0.0075 

12 

7 

8 

-0.12D-02 

2.906 

0.0176 

0.2388 

0.495D-01 

0. 144D-01 -0. 131D+00 

.03422 .01456 

0 

0 .125 

66 .080 

-4 

0.0017 

2 

8 

10 

0.11D-02 

2.848 

0.0178 

0.1992 

0.554D-01 

0. 127D-01 -0. 107D+00 

.03448 .01445 

0 

0 .122 

66 

.040 

-4 

0.0010 

2 

9 

12 

0.99D-03 

2.818 

0.0179 

0.1681 

0.608D-01 

0.120D-01 -0.905D-01 

.03478 .01431 

0 

0 .098 

66 

.032 

-4 

-.0015 

10 

10 

14 

0.920*03 

2.790 

0.0178 

0.1490 

0.647D-01 

0, 119D-01 -0.793D-01 

.03503 .01410 

0 

0 .069 

66 .019 

-4 

-.0019 

10 

11 

16 

0.82D-03 

2.766 

0.0177 

0.1386 

0.674D-01 

0. 123D-01 -0.711D-01 

.03520 .01392 

0 

0 .044 

66 .014 

-16 

-.0021 

10 

12 

18 

0.70D-03 

2.747 

0.0176 

0.1344 

0.693D-01 

0. 129D-01 -0.653D-01 

.03527 .01386 

0 

0 .034 

69 .014 

-16 

-.0021 

10 

13 

20 

0.59D-03 

2.734 

0.0176 

0.1347 

0.707D-01 

0. 137D-01 — 0.617D-01 

.03519 .01383 

0 

0 .030 

69 .016 

-16 

-.0020 

10 

14 

22 

0.49D-03 

2.724 

0.0175 

0.1373 

0.715D-01 

0. 145D-01 -0.594D-01 

.03507 .01380 

0 

0 .029 

65 

.016 

-16 

-.0018 

10 

15 

24 

0.38D-03 

2.719 

0.0174 

0.1411 

0.720D-01 

0. 152D-01 -O.583D-01 

.03492 .01378 

0 

0 .031 

65 

.015 

-16 

-.0015 

10 

16 

26 

0.27D-03 

2.716 

0.0174 

0.1453 

0.724D-01 

0. 1580-01 -0.579D-01 

.03476 .01376 

0 

0 .031 

65 .013 

-16 

-.0013 

9 

17 

28 

0.18D-03 

2.714 

0.0173 

0.1492 

0.725D-01 

0. 163D-01 -0.58CD-01 

.03461 .01374 

0 

0 .030 

65 

.011 

-16 

-.0010 

8 

18 

30 

0.12D-03 

2.713 

0.0173 

0.1528 

0.726D-01 

0. 168D-01 -0.583D-01 

.03448 .01374 

0 

0 .026 

65 

.009 

-16 

-.0008 

7 

19 

32 

0.81D-04 

2.712 

0.0173 

0.1558 

0.726D-01 

0. 172D-01 -0.588D-01 

.03436 .01374 

0 

0 .026 

65 .007 

-16 

-.0007 

6 

20 

34 

0.65D-04 

2.712 

0.0172 

0.1583 

0.726D-01 

0. 175D-01 -0.593D-01 

.03426 .01373 

0 

0 .024 

65 .007 

-3 


6 

21 

36 

0.53D-04 

2.711 

0.0172 

0.1603 

0.725D-01 

0. 177D-01 -0.598D-01 

.03417 .01373 

0 

0 .022 

65 .006 

-3 

SMI 

5 

22 

38 

-0.50D-04 

2.711 

0.0172 

0.1619 

0.725D-01 

0. 179D-01 -0.602D-01 

.03410 .01372 

0 

0 .020 

65 .006 

-3 

-.0004 

5 

23 

40 

-0.51D-04 

2.710 

0.0172 

0.1631 

0.724D-01 

0. 1 BOD-01 -0. 606D-01 

.03404 .01371 

0 

0 .018 

65 .006 

-3 

-.0003 

5 

24 

42 

-C.50D-04 

2.710 

0.0172 

0.1641 

0.724D-01 

0. 181D-01 -0.609D— 01 

.03399 .01371 

0 

0 .016 

65 

.005 

-3 

-.0002 

5 

25 

44 

-0.49D-04 

2.709 

0.0172 

0.1649 

0.723D-01 

0. 182D-01 -0.612D-01 

.03394 .01370 

0 

0 .015 

65 

.005 

-3 

-.0002 

5 

26 

46 

0.20D-04 

2.708 

0.0171 

0.1656 

0.722D-01 

0. 182D-01 -0.614D-01 

.03391 .01370 

0 

0 .013 

65 .005 

-3 

-.0002 

5 


MESH NUHBER 3 RESULTS 


LOWER SURFACE TRANSITION POINT HAS MOVED FROM I- 69 TO I- 
SOURCE STRENGTH HAS BEEN SMOOTHED ISTAG- 79.5 ILE- 82 


70 ITBL- 


70 0.5140D-01 0.55?7D-04 
69 0.6017D-01 0.5929D-04 
68 0, 6958D-01 0.1106D-03 
67 0.7961D-01 0. 1619D-03 


I X OSMTH 

70 0.5140D-01 0. 55770-04 
69 0.6017D-01 0.8318D-04 
68 0.6958D-01 0.1106D-03 
67 0.7961D-01 0. 1619D-03 


X 

CF 

D-STAR 

D-THETA 

H 

HBAR 

Q 

CE 

UEI 

UEV 

0.2010D-02 

C.1000D+11 

0.2043D-04 

0.9340D-05 

2.19 

0.0 

0.8899D-04 

0.0 

0.0 

0.0 

0. 1361D-02 

0.2939D+00 

0. 1128D-04 

0.5542D-05 

2.04 

0.0 

0.8899D-04 

0.0 

0.0785 

0.0 

0.3183D-03 

0.2371D-01 

0.3152D-04 

0. 1492D-04 

2.11 

0.0 

0.8899D-04 

0.0 

0.3493 

0.0 

0.87B6D-05 

0. 1397D-01 

0.2753D-04 

0. 1152D-04 

2.39 

0.0 

0.B899D-04 

0.0 

0.6492 

0.0 

0.4305D-03 

0 . 9461D-02 

0.3003D-04 

0. 1134D-04 

2.65 

0.0 

0. 8899D-04 

0.0 

0.9278 

0.0 

0. 1586D-02 

0. 6806D-02 

0.3711D-04 

0. 1284D-04 

2.89 

0.0 

0.B899D-04 

0.0 

1.1455 

0.0 

0.348CD-02 

0.5080D-02 

0.4770D-04 

0. 1541D-04 

3.10 

0.0 

0.8899D-04 

0.0 

1.2968 

0.0 

0.6112D-02 

0.3884D-02 

0.6090D-04 

0. 1864D-04 

3.27 

0.0 

0.8899D-04 

0.0 

1.3990 

0.0 

0.9478D-02 

0.3015D-02 

0.7631D-04 

0.2237D-04 

3.41 

0.0 

0.8B99D-04 

0.0 

1.4646 

0.0 

0. 1358D-01 

0.2418D-02 

0.9316D-04 

0.2655D-04 

3.51 

0.0 

0.8899D-04 

0.0 

1.5031 

0.0 

0. 1842D-01 

0.1962D-02 

0. 1086D-03 

0.3094D-04 

3.51 

0.0 

0.8899D-04 

0.0 

1.5264 

0.0 

0 . 2399D-01 

0. 1624D-02 

0.1249D-03 

0.3558D-04 

3.51 

0.0 

0.8899D-04 

0.0 

1.5380 

0.0 

0.3028D-01 

0.5247D-02 

0.8585D-04 

0.4052D— 04 

2.12 

1.47 

0.8B99D-04 

0. 1184D-01 

1.5422 

1.5422 

0.3729D-01 

0.4492D-02 

0. 1237D-03 

0.5790D-04 

2.14 

1.48 

0.1282D-03 

0. 1341D-01 

1.5416 

1.5437 

0.4499D— 01 

0.4213D-02 

0. 1622D-03 

0.7724D-04 

2.10 

1.46 

0. 1684D-03 

0.1324D-01 

1.5394 

1.5371 

0.5336D-01 

0.4094D-02 

0.1928D-03 

0.9278D-04 

2.08 

1.43 

0. 1999D-03 

0. 1209D-01 

1.5371 

1.5418 

0.6239D-01 

0.3892D-02 

0.2314D-03 

0.1118D-03 

2.07 

1.43 

0.2399D-03 

0. 1216D-01 

1.5350 

1.5415 

0.7206D-01 

0. 3727D-02 

0.2724D-03 

0.1326D-03 

2.05 

1.42 

0.2829D-03 

0. 1245D-01 

1.5331 

1.5357 

0.8236D-01 

0.3633D-02 

0.3101D-03 

0.1522D-03 

2.04 

1.41 

0.3224D-03 

0. 1223D-01 

1.5316 

1.5324 

0.9326D-01 

0.3554D-02 

0.3491D-03 

0.1725D-03 

2.02 

1.40 

0.3631D-03 

0. 1187D-01 

1.5307 

1.5314 

0.1047D+00 

0. 3481D-02 

0.3887D-03 

0.1930D-03 

2.01 

1.39 

0.4043D-03 

0. 1159D-01 

1.5303 

1.5308 

0.1168D+00 

0. 3411D-02 

0. 4298D-03 

0.2144D-03 

2.00 

1.39 

0.4471D-03 

0* 1137D-01 

1.5304 

1.5306 

0, 1294D+00 

0.3350D-02 

0.4708D-03 

0.2357D-03 

2.00 

1.38 

0.4897D-03 

0. 1116D-01 

1.5310 

1.5310 

0 . 1425D+00 

0 . 3296D-02 

0.5126D-03 

0.2573D-03 

1.99 

1.37 

0. 53300-03 

0. 1093D-01 

1.5322 

1.5321 

0.1562D+00 

0.3247D-02 

0.5544D-03 

0.2790D-03 

1.99 

1.37 

0.5762D-03 

0. 1072D-01 

1.5337 

1.5336 

0.1703D+00 

0. 3199D-02 

C.5977D-03 

0. 3013D-03 

1.98 

1.37 

0.6209D-03 

0. 1054D-01 

1.5355 

1.5352 

0. 1849D+00 

0.3156D-02 

0.6408D-03 

0.3236D-03 

1.98 

1.36 

0.6653D-03 

0. 1038D-01 

1.5376 

1.5371 

0. 1999D+00 

0.3116D-02 

0.6846D-03 

0.3461D-03 

1.98 

1.36 

0.7103D-03 

0. 1022D-01 

1.5398 

1.5393 

0.2154D+00 

0.3080D-02 

0.7285D-03 

0.3687D-03 

1.98 

1.35 

0.7553D-03 

0. 1007D-01 

1.5423 

1.5418 

0.2312D+00 

0.3046D-02 

0.7729D-03 

0.3915D-03 

1.97 

1.35 

0.B007D-03 

0.9909D-02 

1.5450 

1.5445 

0.2474D+00 

0.3015D-02 

0.8173D-03 

0.4142D-03 

1.97 

1.35 

0.8459D-03 

0.9748D-02 

1.5479 

1.5474 


( 5 ) 


56 


0.2640D+00 

0. 2984D-02 

0.8624D-03 

C.4371D-03 

1.97 

0. 2B09D+00 

0. 29550-02 

0.9077D-03 

0. 46010-03 

1.97 

0 . 29810+00 

0.29290-02 

0.95270-03 

0.48280-03 

1.97 

0. 3156D+00 

0.29050-02 

0. 99770-03 

G.5054D-03 

1.97 

0. 3333D+00 

• 0.2882D-02 

0. 1043D-C2 

0.52780-03 

1.98 

C.3512D+00 

0.2S60D-C2 

C. 1OS7D-02 

0.550OD-03 

1.98 

0.3B94D+0O 

C. 2S38D-02 

C. 1122D-02 

O.5721D-03 

1.98 

0.3876D+00 

0 . 2816D-02 

0. 11780-02 

0.59450-03 

1.98 

0 . 4061D+00 

0.2798D-02 

0. 1222D-02 

0.6158D-03 

1.98 

0 . 4246D+00 

0.2782D-02 

0. 1264D-02 

0.6359D-03 

1.99 

0.4432D+00 

0.2764D-02 

0. 1309D-02 

0.6574D-03 

1.99 

0 . 4619D+00 

0.2745D-02 

0. 1351D-02 

0.6767D-03 

2.00 

0 . 4B06D+00 

0.2738D-02 

0. 1399D-02 

0.6989D-03 

2.00 

0. 4993D+00 

0.2731D-O2 

0. 1416D-02 

0.7049D-03 

2.01 

0.5179D+00 

0.2685D-02 

0. 1502D-02 

0. 74580-03 

2.01 

0.53660+00 

0.26000-02 

0. 1528D-02 

0.75700-03 

2.02 

0.5551D+00 

0.2304D-C2 

0. 1784D-02 

0.88020-03 

2.03 

0.5736D+00 

0. 1433D-02 

0.22220-02 

0. 1009D-02 

2.20 

0.5921D+00 

0.5418D-03 

0.4393D-02 

0. 1641D-02 

2.68 

0.6104D+00 

0.33390-03 

0.5317D-02 

0.18160-02 

2.93 

0.6286D+00 

0. 3678D-03 

0. 53020-02 

0. 1866D-02 

2.84 

0.6466D+00 

0. 43900-03 

0.56160-02 

0.20810-02 

2.70 

0.6645D+00 

0.53750-03 

0.54970-02 

0.21520-02 

2.55 

0. 68230+00 

0.6326D-03 

0.55950-02 

0.2302D-02 

2.43 

0. 69980+00 

0.7156D-03 

0. 56250-02 

0.24070-02 

2.34 

0.7170D+00 

0.7874D-03 

0.5697D-02 

0.2519D-02 

2.26 

0. 73400+00 

0.84220-03 

0.58040-02 

0.2631D-02 

2.21 

0.7507D+00 

0.88160-03 

0.59300-02 

0.27400-02 

2.16 

0.76700+00 

0.9052D-03 

0.61120-02 

0.2B64D-02 

2.13 

C.7830D+00 

0.91760-03 

0.62950-02 

0.29790-02 

2.11 

0.7987D+00 

0.9177D-03 

0.6538D-02 

0.31130-02 

2.10 

0.8139D+00 

0.9089D-03 

0.67860-02 

0.3240D-02 

2.09 

0.82B8D+00 

0.8913D-03 

0.7089D-02 

0.33850-02 

2.09 

0.8431D+00 

0.8677D-03 

0.7395D-02 

0.35230-02 

2.10 

0.857CD+00 

0.8379D-03 

0.7757D-02 

0.3678D-02 

2.11 

0.8704D+00 

0.80410-03 

0.8121D-02 

0.38250-02 

2.12 

0.8833D+00 

0.7656D-03 

0.8543D-02 

0.3988D-02 

2.14 

0.89560+00 

0.7249D-03 

0.8966D-02 

0.4142D-02 

2.16 

0.9074D+00 

0.6814D-03 

0.9444D-02 

0.4309D-02 

2.19 

0.91650+00 

0.63700-03 

0.99250-02 

0.44650-02 

2.22 

0.9291D+00 

0.59190-03 

0. 1045D-01 

0.4628D-02 

2.26 

0.9390D+00 

0.5473D-03 

0. 1097D-01 

0.47810-02 

2.30 

0.94820+00 

0.50360-03 

0.1152D-01 

0.49320-02 

2.34 

0.9568D+00 

0.4617D-03 

0. 1207D-01 

0. 50720-02 

2.38 

0.96460+00 

0.4220D-03 

0.1262D-01 

0.5203D-02 

2.43 

0.9717D+00 

0.38480-03 

0. 1317D-01 

0.53250-02 

2.47 

0.9781D+00 

0. 35040-03 

0.13690-01 

0.5431D-02 

2.52 

0.9837D+00 

0.3186D-03 

0.14210-01 

0.5 530D-02 

2.57 

0.9885D+00 

0.28970-03 

0. 1469D-01 

0.56090-02 

2.62 

0.9925D+00 

0 . 2629D-03 

0. 1516D-01 

0.56820-02 

2.67 

0.9957D+00 

0.2379D-03 

0. 1558D-01 

0.57350-02 

2.72 

0.9981D+00 

0.21360-03 

0. 1601D-01 

0.577 9 D-02 

2.77 

0.9995D+00 

0. 1891D-03 

0. 1643D-01 

0.58C9D-02 

2.83 

0.1000D+01 

0. 1644D-03 

0. 1686D-01 

0.58270-02 

2.89 

0. 1001D+01 

0.0 

0. 1683D-01 

C.5B39D-02 

2.88 

0. I002D+01 

0.0 

0. 1673D-01 

0.5848D-02 

2.86 

0.1003D+01 

0.0 

0.16580-01 

0.5B58D-02 

2.83 

0. 1006D+01 

0.0 

0. 1640D-01 

0.5865D-02 

2.79 

0. 1009D+01 

0.0 

0.16130-01 

0.59000-02 

2.73 

0.1014D+01 

0.0 

0.15850-01 

0.5934D-02 

2.67 

0. 1019D+01 

0.0 

0.15450-01 

0.5949D-02 

2.60 

0. 1026D+01 

0.0 

0.1503D-01 

0.59820-02 

2.51 

0. 1034D+01 

0.0 

0.14470-01 

0.5978D-02 

2.42 

0. 1043D+01 

0.0 

0. 1389D-01 

0.5980D-02 

2.32 

0. 1054D+C1 

0.0 

0. 1322D-01 

0.5949D-02 

2.22 

0.1067D+01 

0.0 

0.12520-01 

0.5897D-02 

2.12 

0.1082D+01 

0.0 

0.1183D-01 

0.5833D-02 

2.03 

0.U00D+01 

0.0 

0.11120-01 

0.5735D-02 

1.94 

0. 11210+01 

0.0 

0.10460-01 

0.5635D-02 

1.86 

0. 1146D+01 

0.0 

0.9792D-02 

0.5499D-02 

1.78 

0. 1175D+01 

0.0 

0.9260D -02 

0.54010-02 

1.71 

0. 1209D+01 

0.0 

0.87340-02 

0.52760-02 

1.66 

0.12500+01 

0.0 

0.8288D-02 

0.5173D-02 

1.60 

0.12980+01 

0.0 

0.78430-02 

0.50450-02 

1.55 

0. 1357D+01 

0.0 

0.75190-02 

0.4973D-02 

1.51 

0. 1430D+01 

0.0 

0.7185D-02 

0.4875D-02 

1.47 

0.15200+01 

0.0 

0.6908D-02 

0.4802D-02 

1.44 

0. 16350+01 

0.0 

0. 6647D-02 

0.4725D-02 

1.41 

0. 1785D+01 

0.0 

0.64240-02 

0.46630-02 

1.38 

0. 1988D+01 

0.0 

0.62060-02 

0.4595D-02 

1.35 

0. 22760+01 

0.0 

0.6032D-02 

0.45510-02 

1.33 

0. 2716D+01 

0.0 

0.5876D-02 

0.4513D-02 

1.30 

0.3465D+01 

0.0 

0.5721D-02 

0.4472D-02 

1.28 

0.5032D+01 

0.0 

0.5561D-02 

0.4424D-02 

1.26 

0.1047D+02 

0.0 

0.5379D-02 

0.43730-02 

1.23 

X-IMAX- 0.574 IMAX-127 

DUMAX- 0.008 

DUAVG- 0.002 

1STAG* 79.5 

ILE- B2 





1.34 

0. 8916D-03 

0.9603D-02 

1.5508 

1.5504 

1.34 

0.93750-03 

0.9463D-02 

1.5539 

1.5535 

1.34 

0.98290-03 

0.9306 D-02 

1.5572 

1.5570 

1.34 

0. 1C28D-C2 

0.9125D-02 

1.5609 

1.5607 

1.33 

0. 1073D-C2 

0. 8927D-C2 

1.5647 

1.5647 

1.33 

0. 11170-02 

0.E73CD-02 

1.5688 

1.5690 

1.33 

C. 1162D-02 

C.8559D-C2 

1.5729 

1.5732 

1.33 

0.12070-02 

0.8388D-02 

1.5774 

1.5774 

1.33 

0. 12500-02 

0. 81540-02 

1.5819 

1.5824 

1.32 

0. 12900-02 

0.7918D-02 

1.5868 

1.5878 

1.32 

0.13330-02 

0.76750-02 

1.5923 

1.5932 

1.32 

0.13730-02 

0.74280-02 

1.5990 

1.5981 

1.32 

0. 1416D-02 

0.67850-02 

1.6056 

1.6072 

1.32 

0.14300-02 

0.65650-02 

1.6119 

1.6136 

1.32 

0. 1515D-02 

0.6679D-02 

1.6149 

1.6172 

1.34 

0. 1551D-02 

0.87910-02 

1.6073 

1.6016 

1.39 

0.18460-02 

0. 1426D-01 

1.5443 

1.5478 

1.66 

0.2377D-02 

0.24920-01 

1.4029 

1.3954 

2.14 

0.4717D-02 

0.3127D-01 

1.2873 

1.2914 

2.39 

0.5695D-02 

0.32B 10-01 

1.2463 

1.2523 

2.34 

0.56640-02 

0.3384D-01 

1.2220 

1.2263 

2.23 

0.5980D-02 

0.34620-01 

1.2026 

1.2039 

2.12 

0.58350-02 

0.34940-01 

1.1843 

1.1849 

2.02 

0.59200-02 

0.35000-01 

1.1671 

1.1685 

1.95 

0.5932D-02 

0.34810-01 

1.1549 

1.1535 

1.89 

0.59890-02 

0.3429D-01 

1.1410 

1.1399 

1.85 

0. 60800-02 

0.33790-01 

1.1283 

1.1269 

1.82 

0.6190D-02 

0.33390-01 

1.1158 

1.1143 

1.80 

0.63570-02 

0.3307D-01 

1.1037 

1.1018 

1.79 

0.6522D-02 

0.3281D-01 

1.0914 

1.0695 

1.79 

0.67460-02 

0.32650-01 

1.0791 

1.0772 

1.79 

0.6971D-02 

0.3259D-01 

1.0671 

1.0649 

1.80 

0.7250D-02 

0.3263D-01 

1.0551 

1.0528 

1.81 

0.75280-02 

0.32770-01 

1.0433 

1.0409 

1.82 

0.7858D-02 

0.3299D-01 

1.0318 

1.0293 

1.84 

0.81870-02 

0.3328D-01 

1.0206 

1.0179 

1.87 

0.85700-02 

0.33620-01 

1.0098 

1.0070 

1.89 

0. 89510-02 

0.3401D-01 

0.9993 

0.9964 

1.92 

0.93820-02 

0.3444D-01 

0.9894 

0.9863 

1.96 

0.98140-02 

0.34890-01 

0.9800 

0.9768 

1.99 

0.10280-01 

0.35360-01 

0.9712 

0.9679 

2.03 

0. 1076D-01 

0.35610-01 

0.9630 

0.9597 

2.07 

0.11250-01 

0.35810-01 

0.9556 

0.9522 

2.12 

0.11740-01 

0.36000-01 

0.948B 

0.9454 

2.16 

0.12230-01 

0.3617D-01 

0.9427 

0.9394 

2.21 

0. 1272D-01 

0.3633D-01 

0.9374 

0.9341 

2.26 

0.1319D-01 

0.3647D-01 

0.9327 

0.9296 

2.31 

0.13660-01 

0.3659D-01 

0.9287 

0.9256 

2.35 

0. 1409D— 01 

0.3666D-01 

0.9253 

0.9224 

2.40 

0. 1452D-01 

0.36680-01 

0.9225 

0.9197 

2.45 

0. 1491D-01 

0.3667D-01 

0.9203 

0.9176 

2.50 

0. 1530D-01 

0.3666D-01 

0.9185 

0.9160 

2.55 

0. 1570D-01 

0.36660-01 

0.9168 

0.9149 

2.61 

0.16100-01 

0.36680-01 

0.9180 

0.9143 

2.60 

0. 1606D-01 

0.3667D-01 

0.9178 

0.9142 

2.58 

0. 1597D-01 

0.3665D-01 

0.9172 

0.9139 

2.55 

0.1582D-01 

0.3661D-01 

0.9164 

0.9134 

2.51 

0.15650-01 

0.3657D-01 

0.9153 

0.9127 

2.47 

0. 1538D-01 

0.3655D-01 

0.9139 

0.9119 

2.41 

0.15100-01 

0.3656D-01 

0.9125 

0.9110 

2.34 

0. 14710-01 

0.3660D-01 

0.9109 

0.9101 

2.26 

0.1431D-01 

0.366BD-01 

0.9098 

0.9094 

2.18 

0. 1377D-01 

0.36820-01 

0.9088 

0.9090 

2.09 

0. 1322D-01 

0.3702D-01 

0.9085 

0.9092 

1.99 

0. 1259D-01 

0.37320-01 

0.9090 

0.9102 

1.90 

0.1194D-01 

0.37790-01 

0.9104 

0.9121 

1.81 

0. 11300-01 

0.38220-01 

0.9132 

0.9151 

1.72 

0.10650-01 

0.3857D-01 

0.9172 

0.9190 

1.64 

0. 1004D-01 

0.3B82D-01 

0.9220 

0.9241 

1.57 

0.94370-02 

0.38960-01 

0.9274 

0.9299 

1.51 

0.8956D-02 

0.3898D-01 

0.9337 

0.9359 

1.45 

0.8478D-02 

0 .38840-01 

0.9402 

0.9422 

1.40 

0.8074D-02 

0.38480-01 

0.9465 

0.9488 

1.35 

0.76680-02 

0.3792D-01 

0.9527 

0.9553 

1.31 

0.7375D-02 

0.37160-01 

0.9590 

0.9612 

1.27 

0.706BD-02 

0.36140-01 

0.9649 

0.9668 

1.24 

0.68150-02 

0.3484D-01 

0.9705 

0.9722 

1.21 

0.65740-02 

0.3325D-01 

0.9757 

0.9773 

1.18 

0.63690-02 

0.3135D-01 

0.9805 

0.9822 

1.15 

0.61670-02 

0.2913D-01 

0.9850 

0.9868 

1.13 

0.6005D-02 

0.2659D-01 

0.9889 

0.9906 

1.11 

0.58590-02 

0.23670-01 

0.9924 

0.9937 

1.08 

0.5713D-02 

0.2023D-01 

0.9953 

0.9971 

1.06 

0.5566D-02 

0.16030-01 

0.9976 

1.0017 

1.04 

0.53750-02 

0. 1081D-01 

0.9992 

0.9988 


LOWER LIMIT ON LAMBDA EXCEEDED ON UPPER SURFACE 


I 

K 

PIW 

DELTA 

163 

1 

0.300000D+01 

0.0 

163 

2 

0.1655150+02 

0.1701360+00 

163 

3 

0.1934490+02 

0 . 348077D-01 

163 

4 

0. 193648D+02 

0.3345030-01 

163 

5 

0. 193652D+02 

0.334177D-01 


USER'S LOWER SURFACE TRANSITION POINT TOO CLOSE TO STAGNATION POINT 
RESET TO XTL-0.036 FOR THIS MESH 



USER’S LOWER SURFACE TRANSITION POINT TOO CLOSE TO STAGNATION POINT 
RESET TO XTL-0.043 FOR THIS MESH 


USER'S LOHER SURFACE TRANSITION POINT TOO CLOSE TO STAGNATION POINT 
RESET TO XTL-0.051 FOR THIS MESH 


LOHER SURFACE 

0.2010D-02 
0.3145D-02 
0.5665D-02 
0.89120*02 
0.1288D-01 
0. 1757D-01 
0. 22960-01 
0. 2905D-01 
0.3583D-01 
0.43280-01 
0.5140D-01 
0.60170-01 
0 . 6958D-01 
0.7961D-01 
0. 90240-01 
0. 1G15D+00 
0. 1132D+00 
0 . 1256D+C0 
0. 1384D+00 
0. 1518D+00 
0 . 1656D+00 
0. 1799D+00 
0. 1946D+00 
0. 2097D+00 
C.2251D+00 
0. 2410D+00 
0.2571D+00 
0.2735D+C0 
G.2901D+C0 
0.3070D+00 
0.3240D+G0 
0 . 3412D+G0 
0.35860+00 
0.3761D+00 
0. 39370+00 
0.4115D+00 
0.4295D+00 
0. 4476D+00 
0.4659D+00 
0. 4843D+00 
0.50290+00 
0.5216D+00 
C.5405D+00 
0.5594D+00 
0.57840+00 


CF 

0 . 10000+11 
0.11840+00 
0. 15950-01 
0.99950-02 
0. 7126D-02 
0.54480*02 
0.43620-02 
0. 36080-02 
0, 3059D-02 
0.2649D-02 
0.8004D-02 
0.5114D-02 
0. 4886D-02 
0.4639D-G2 
G.4617D-02 
0 . 4419D-02 
0.42140-02 
0.41120-02 
0. 40230-02 
0.39360-02 
0.38640-02 
0.3795D-02 
0.37260-02 
0.3666D-02 
0.3612D-02 
0.3566D-02 
0.3527D-02 
0.3492D-02 
0. 3461D-02 
0.3428D-02 
0.3385D-02 
0.33320-02 
0.3263D-02 
0. 31860-02 
0.3089D-02 
0.2991D-02 
0.2887D-02 
0. 2788D-02 
0.2683D-02 
0.25900-02 
0.25030-02 
0.2426D-02 
0.23580-02 
0.2300D-02 
0.22510-02 


D-STAR 
0 . 2043D-04 
0.16C5D-04 
0.4588D-04 
0.4968D-04 
0.56950*04 
0.6565D-04 
0.7516D-04 
0.85270-04 
0.9476D-04 
0. 1045D-03 
0.6327D-04 
0.11440-03 
0. 1327D-03 
0, 1676D-03 
0. 1821D-03 
0. 2108D-03 
0.24400-03 
0.2709D-03 
0.2991D-03 
0.3287D-03 
0.3571D-03 
0.3862D-03 
0. 4168D-03 
0.44680-03 
0.4770D-G3 
0.5068D-03 
0.53490-03 
0.56290*03 
0.58920-03 
0.6148D-03 
0.64340-03 
0.67290-03 
0.71000-03 
0.74750-03 
0.7997D-03 
0.B532D-03 
0.9194D-03 
0.9873D-03 
0. 1073D-02 
0.1159D-02 
0.12570-02 
0.1356D-02 
0.1464D-02 
0. 1573D-02 
0. 16B6D-02 


D-THETA 
0. 9340D-05 
0.78730-05 
0. 1977D-04 
0.20370-04 
0.22810-04 
0.2583D-04 
0.2915D-04 
0.32650-04 
0.36280-04 
0.4001D-04 
0. 43770-04 
0.7041D-04 
0.8184D-04 
0. 1045D-03 
0.1143D-03 
0. 1322D-03 
0. 1525D-03 
0.16960-03 
0.1875D-03 
0.2060D-03 
0.22380-03 
0.24190-03 
0.2609D-03 
0.27940-03 
0.2979D-03 
0.31620-03 
0.33330-03 
0. 35020-03 
0.3658D-03 
0.3806D-03 
0.39680-03 
0.4129D-03 
0.43320-03 
0.4534D-03 
0.4820D-03 
0.51140-03 
0.5485D-03 
0.58670-03 
0. 63540-03 
0.6851D-03 
0.74170-03 
0.8004D-03 
0.8647D-03 
0.93050-03 
0.9990D-03 


0.56090-04 
0. 56C9D-04 
0.5609D-04 
0.5609D-04 
0.56090-04 
0.56090-04 
0.5609D-04 
0.56090-04 
0.5609D-04 
0.56090-04 
0.5609D-04 
0.10390-03 
0.12200*03 
0. 1564D-03 
0.17270*03 
0.20280-03 
0. 2367D-03 
0. 265CD-03 
0.2949D-03 
0. 32630-03 
0.3568D-03 
0.38820-03 
0.4212D-03 
C.4538D-03 
0.4867D-03 
0.5195D-O3 
0.55070-03 
0.58210-03 
0.61200-03 
0.6413D-03 
0.6736D-03 
0.70660-03 
0.7472D-03 
0.7876D-03 
0.84270-03 
0. 89840-03 
0.96660-03 
0.1035D-02 
0.11210-02 
0. 1206 D-02 
0.13010-02 
0. 1396D-02 
0.14980-02 
0. 1599D-02 
0. 1702D-02 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

-0.9000D-02 
0. 10760-01 
0. 12140-01 
0. 1123D-01 
0.1019D-01 
0.9975D-02 
0.10320-01 
0. 1021D-01 
0. 99460-02 
0. 9749D-02 
0.9538D-02 
0.9358D-02 
0. 92260-02 
0. 90850-02 
0.8921D-02 
0.87100-02 
0.8429D-02 
0 ,80840-02 
0.7659D-02 
0.721CD-02 
0.6824D-02 
0.66000-02 
0.6600D-02 
0.6885D-02 
0.7477D-02 
0.8280D-02 
0. 9218D-02 
0. 1027D-01 
0. 11400-01 
0. 1250D-01 
0. 13540-01 
0. 1452D-01 
0. 1542D-01 
0. 1621D-01 
0. 1691D-01 


0.5974D+00 0.2207D-02 0.1800D-02 0.1070D-02 
0.6164D+00 0.2169D-02 0.1919D-02 0.1143D-02 
0.63550+00 0.21340-02 0.2039D-02 0.1218D-02 
0. 6544D+00 0.2101D-02 0.2164D-02 0.1297D-02 
0. 6733D+00 0.2071D-02 0.2291D-C2 0.1376D-02 
0. 69200+00 0. 20430-02 0.2421D-C2 0.1458D-02 
0.7106D+00 0. 20190-02 0.25510-02 0.1540D-02 
0.7290D+00 0.19960*02 0.26800*02 0.16230-02 
0.74720+00 0. 1978D-02 0.2811D-02 0.1707D-02 
0.7651D+00 0, 1965D-02 0.2929D-02 0. 17850-02 
0.7826D+00 0.1956D-02 0.3047D-02 0.1863D-02 
0.799BD+00 0. 1950D-02 0.31550-02 C.1936D-02 
0.B165D+CD 0.1946D-B2 0.3264D-C2 0.2C10O-02 
0. 8328D+00 0.1945D-02 C.3360D-02 0.2077D-02 
0.8485D+00 0. 1946D-02 0.34550-02 0.2144D-02 
0. 8637D+00 C.1952D-02 0.3531D-02 C.2201D-C2 
0.8784D+00 0. 1962D-02 0.3605D-02 0.2258D-02 
0.8923D+00 0. 1978D-02 0.3650D-02 0.2298D-C2 
0. 90560+00 0. 1999D-02 0.3691D-02 0.2337D-02 
0.9162D+00 0.2028D-02 0.3697D-02 0.23560-02 
0.9299D+00 0.20620*02 0.36970-02 0.2373D-C2 
C.9409D+00 0.21C5D-02 0.3661D-02 0.2367D-02 
0. 95100+00 0.2153D-02 0.3620D-02 0.23600-02 
0.9602D+00 0.221CD-02 0. 35420-02 0.2329D-02 
0.96840+00 C.2272D-02 0.3457D-02 0.2293D-02 
0.9758D+00 0.23410-02 0.3349D-02 0.224OD-C2 
0.9822D+00 0.24100-02 0.32400-02 0.2186D-02 
0.9B76D+00 0.2477D-02 0.31310-02 0.2128D-02 
0.99200+00 0. 25410-02 0.3027D-02 0.2071D-02 
0.9955D+00 0.25990-02 0.29310-02 0.2C17D-02 
0.9980D+00 0.26450-02 0.2860D-02 0.1976D-02 
0.9995D+C0 0.26900-02 0.279 0D-02 0.1935D-02 
0. 1000D+C1 0.2741D-02 0.2699D-02 0.1879D-C2 
0. 1001D+01 0.0 0.2682D-02 0.18670-02 
0. 1002D+01 0.0 0.26BBD-02 0.1872D-C2 
0. 1003D+01 0.0 0.2701D-02 0.1881D-02 
0. 1006D+01 0.0 0.2708D-02 0.1886D-02 
0.10090+01 0.0 0.2725D-02 0.18990-02 
0. 1014D+01 0.0 0.2742D-02 0.1912D-02 
0. 1019D+01 0.0 0.2748D-02 0. 19180-02 
0.1026D+01 0.0 0.2750D-02 0.1922D-02 
0. 1034D+01 0.0 0,27590-02 0.1933D-02 
0.1043D+01 0.0 0.2761D-02 0.1938D-02 
0.1054D+01 0.0 0.27450-02 0.1934D-02 
0. 1067D+01 0.0 0.2733D-02 0.1933D-C2 
0.10820+01 0.0 0.2679D-02 0.1906D-02 
0. 1100D+01 0.0 0.2624D-02 0.1878D-02 
0. 11210+01 0.0 0.2565D-02 0.1848D-02 
0. 1146D+01 0.0 0.2502D-02 0.1816D-02 
C. 1175D+01 0.0 0.2428D-02 0.1775D-02 
0.12090+01 0.0 0.2360D-02 0.1739D-02 
0.12500+01 0.0 0.2297D-02 0.1706D-02 
0. 1298D+01 0.0 0.2233D-02 0.1672D-02 
0. 1357D+01 0.0 0.2175D-02 0.16410-02 
0.14300+01 0.0 0.2117D-02 0.1610D-02 
0.15200+01 0.0 0.2068D-02 0.15B6D-02 
0.1635D+01 0.0 0.20200-02 C.1562D-02 
0.17850+01 0.0 0.1975D-02 0.1540D-02 
0. 1988D+01 0.0 0. 1932D-02 0.15190-02 
0. 22760+01 0.0 0.18950-02 0.1503D-C2 


1.68 

1.45 

0.1805D-02 

0.1754D-01 

1.0088 

1.0059 

1.68 

1.45 

0.19100-02 

0.1811D-01 

0.9938 

0.9909 

1.67 

1.45 

0.2016D-02 

0. 1864D-01 

0.9794 

0.9765 

1.67 

1.45 

0.2124D-02 

0. 1915D-01 

0.9654 

0.9625 

1.66 

1.46 

0.2232D-02 

0.1963D— 01 

0.9519 

0.9490 

1.66 

1.46 

0.2341D -02 

0.2008D-01 

0.9390 

0.9361 

1.66 

1.46 

0. 2449D-02 

0. 20490-01 

0.9267 

0.9239 

1.65 

1.46 

0.2556D-02 

0.2089D-01 

0.9152 

0.9123 

1.65 

1.46 

0.26620-02 

0.21250-01 

0.9044 

0.9015 

1.64 

1.46 

0.2757D-02 

0. 21570-01 

0.8944 

0.8916 

1.64 

1.46 

0.2851D-C2 

0.2187D-01 

0.8853 

0.8826 

1.63 

1.46 

0.2935D-02 

G.2213D-01 

0.8771 

0.8744 

1.62 

1.45 

0.3021D-02 

0.22380-01 

0.8694 

0.6668 

1.62 

1.45 

0.3094D-C2 

0.2262D-01 

0.8624 

0.8598 

1.61 

1.45 

0. 31670-02 

0.2285D-01 

0.8561 

0.8537 

1.60 

1.44 

0.3224D-02 

0. 23060-01 

0.8507 

0.8483 

1.60 

1.44 

0.3281D-02 

0.2325D-01 

0.8462 

0.8438 

1.59 

1.43 

0.3313D-C2 

0.2340D-01 

0.8426 

0.8403 

1.58 

1.42 

0.3344D-02 

0.2352D-01 

0.8400 

0.8380 

1.57 

1.41 

0.3347D-02 

0.23590-01 

0.8387 

0.8367 

1.56 

1.40 

0.3347D-C2 

0.23590-01 

0.8385 

0.8367 

1.55 

1.39 

0.3317D-C2 

0.2353D-01 

0.8396 

0.8380 

1.53 

1.38 

0.3287D-C2 

0.2340D-01 

0.8420 

0.8407 

1.52 

1.37 

0.3225D-02 

0.2317D-01 

0.8459 

0.8448 

1.51 

1.35 

0.3161D-02 

0.2281D-01 

0.8513 

0.8504 

1.49 

1.34 

0.3078D-02 

0.22340-01 

0.8580 

0.8573 

1.48 

1.32 

0.2995D-02 

0.2179D-01 

0.8655 

0.8652 

1.47 

1.31 

0.2910D-02 

0.2112D-C1 

0.8736 

0.8735 

1.46 

1.30 

0.2831D-02 

0.2033D-01 

0.BB17 

0.8819 

1.45 

1.29 

0.2756D-02 

0. 1932D-01 

0.6896 

0.8899 

1.45 

1.28 

0.2701D-02 

0. 1788D-01 

0.8965 

0.8969 

1.44 

1.27 

0.26470-02 

0.1676D-01 

0.9034 

0.9041 

1.44 

1.26 

0.2574D-02 

0. 1610D-01 

0.9180 

0.9119 

1.44 

1.26 

0.2558D-02 

0. 1609D-01 

0.9178 

0.9125 

1.44 

1.26 

0.25640-02 

0. 1612D-01 

0.9172 

0.9122 

1.44 

1.26 

0.25740-02 

0. 1620D-01 

0.9164 

0.9111 

1.44 

1.26 

0.25790-02 

0. 1630D-01 

0.9153 

0.9098 

1.44 

1.26 

0. 25920-02 

0. 1643D-01 

0.9139 

0.9081 

1.43 

1.26 

0.2605D-02 

0. 1659D-01 

0.9125 

0.9063 

1.43 

1.26 

0.26090-02 

0.1676D-01 

0.9109 

0.9052 

1.43 

1.26 

0.26090-02 

0. 1694D-01 

0.9098 

0.9042 

1.43 

1.26 

0.2616D-02 

0.1716D-01 

0.9088 

0.9030 

1.42 

1.25 

0.26 160-02 

0.17390-01 

0.9085 

0.9023 

1.42 

1.25 

0.2602D-02 

0.1762D-01 

0.9090 

0.9024 

1.41 

1.24 

C. 25920-02 

0.17850-01 

0.9104 

0.9037 

1.41 

1.23 

0.2546D-G2 

0. 18040-01 

0.9132 

0.9070 

1.40 

1.22 

0.25010-02 

0. 1819D-01 

0.9172 

0.9116 

1.39 

1.21 

0.24520-02 

0.18330-01 

0.9220 

0.9166 

1.38 

1.20 

0.2400D-C2 

0 . 1844D-01 

0.9274 

0.9223 

1.37 

1.19 

0.2339D-02 

0. 184BD-01 

0.9337 

0.9289 

1.36 

1.18 

0.22S2D-02 

0. 18480-01 

0.9402 

0.9356 

1.35 

1.17 

0.22290-02 

0. 1842D-01 

0.9465 

0.9421 

1.34 

1.15 

0.21750-02 

0. 1828D-01 

0.9527 

0.9465 

1.33 

1.14 

0.2126D-02 

0.1805D-01 

0.9590 

0.9548 

1.31 

1.13 

0. 2076D-02 

0. 1771D-01 

0.9649 

0.9609 

1.30 

1.12 

0.20340-02 

0. 1725D-01 

0.9705 

0.9665 

1.29 

1.11 

0.19920-02 

0.16650-01 

0.9757 

0.9716 

1.28 

1.09 

0.1953D-02 

0, 1588D-01 

0.9805 

0.9767 

1.27 

1.08 

0.19140-02 

0. 1493D-01 

0.9850 

0.9813 

1.26 

1.07 

0. 1882D-02 

0. 1379D-01 

0.9889 

0.9852 


58 


G.2716D+01 0.0 
0.34650+01 0.0 
0.5032D+01 0.0 
0. 1G47D+02 0.0 


0. 16620-02 0. 14890-02 1.25 1.06 0.1852D-02 0.1242D-01 0 9924 a 9RR#, 

0.18280-02 0,14760-02 1.24 1.05 SillllD-M 0:10750-01 0 9919 

0»1794D-02 0.1462D-02 1.23 1.04 0.1791D-02 0.8658D-02 0.9976 0*9957 

0. 1757D-02 0, 1450D-02 1.21 1.02 0. 17530-02 0.5715D-02 0.9992 o!9943 


1.067 w JMAX* 22 DUMAX- 0,007 DUAVG- 0.003 
I ST AG* 79.5 ILE« 82 


UPPER LIMIT OH LAMBDA EXCEEDED ON LONER SURFACE 


PIH 

0.300000D+01 
0.132903D+01 
0.1058610+01 
0.996900D+00 
0. 981683D+00 
0.977858D+00 
0.976892D+00 
0.976648D+00 


0 0° ELTA 
0I869482D-02 
0. 1149800-01 
0.130084D-01 
0. 134104D-01 
0.135134D-01 
0. 135395D-01 
0.135461D-01 


COIHF- 0. 1165D-01 COSY- 0.1141D-01 

CDINFT- 0.87460-02 CDINFB- 0.2901D-02 

CDF- 0.5950D-02 COFT- 0. 32280-02 CDFB- 0.2723D-02 


ITERATIVE SOLUTION 


ANG DIVISIONS RAD DIVISIONS 
140 32 

CYCLE MAX CORRECN I J A 

1 -0.6286D-04 128 34 

2 0.4213D-04 3 3 

2 0.0 00 

i v S RESIDAL 1 AVG RESIDAL 2 
0.2368638D-03 0. 1783229D-06 


AVG CORRECN MAX RESIDAL 
0.6303D-05 -0.9488D-04 
0.4565D-05 -0.5426D-04 
0.0 0.4235D-04 


WORK 

77.5469 


REDUCTN/HORK 

0.9114 


J AVG RESIDAL CIRC/ALFA HORK 

34 0.3521D-06 0.0848 75.7773 

33 0.2983D-06 0.0848 77.5469 

28 0.1783D-06 0.0848 77.5469 

REDUCTN/CYCLE 

0.8795 


MESH NO.- 3 


BL ITER NO.- 


X- 0.57363 AT THE SHOCK 
CP- -0.84878 UPSTREAM 

CP- -0.48826 DOWNSTREAM, RANKINE-HUGONIOT 
CP- -0.47624 DOWNSTREAM; ISENTROPIC 


ANG OF ATTACK 
2.7063 


MACH NO 
0.7300 


CL(CIRCULATN) 

0.7977 


CL - 0.60308 CD - 0.0166 COB - 0.01661 CM - -0.0875 
COORDINATES OF AIRFOIL HAPPED FROM CIRCLE AND PRESSURE DISTRIBUTION 


ARC LENGTH 
—1,0184 
-1,0179 
-1.0164 
-1,0139 
-1.0104 
-1.0060 
-1.0005 
-0.9942 
-0.9868 
-0.9785 
-0.9693 
-0.9592 
-0.9483 
-0.9365 
-0.9240 
-0.9107 
-0.8966 
-0.8820 


0.0000 

0.0001 

0.0003 

0.0004 

0.0006 

0.0008 

0.0010 

0.0011 

0.0012 

0.0012 

0.0010 

0.0007 

0.0003 

-0,0004 

- 0.0012 

- 0.0022 

-0.0034 


HACH NO 
-0.6862 
-0.6677 
-0.6544 
-0.6460 
-0.6382 
-0,6311 
-0.6243 
-0.6182 
-0.6129 
-0.6086 
-0.6055 
-0.6036 
-0.602B 
-0.6028 
-0.6039 
-0.6058 
-0.6086 
-0.6121 


V/V0 

-0.9454 

-0.9220 

-0.9050 

-0.8943 

-0.8843 

-0.8752 

-0.8664 

-0.8586 

-0.8518 

-0.8463 

-0.8423 

-0.8398 

-0.8387 

-0.8388 

-0.8401 

-0.8426 

-0.8462 

-0,8508 


DENSITY 

1.0286 

1.0404 

1.0489 

1.0542 

1.0591 

1.0635 

1.0678 

1.0715 

1,0747 

1.0774 

1.0792 

1.0804 

1.0809 

1.0809 

1.0803 

1.0791 

1.0774 

1.0752 


HASS FLOW 
-0.9724 
-0.9593 
-0.9493 
-0.9428 
-0.9366 
-0.9308 
-0.9251 
-0.9200 
-0.9155 
-0.9117 
-0,9090 
-0.9073 
-0.9065 
-0.9066 
-0.9075 
-0.9092 
-0.9117 
-0.9148 
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~ 0.8667 

0 . 84 B 5 

- 0.0048 

- 0.8509 

0.8328 

- 0.0064 

- 0.8345 

0.8165 

- 0.0081 

- 0. 8177 

0.7998 

- 0.0100 

- 0.8004 

0.7826 

- 0.0121 

- 0. 7827 

0.7651 

- 0.0143 

- 0.7647 

C . 7^72 

- 0.0167 

- 0.7463 

0.7290 

- 0.0192 

- 0.7278 

0.7106 

- C . C 217 

- 0.7090 

0.6920 

- 0.0244 

- 0.6900 

0.6733 

- 0.0271 

- 0.6710 

0.6544 

- 0.0299 

- 0.6518 

0.6355 

- 0.0326 

- 0.6326 

0.6164 

- 0.0354 

- 0.6134 

0.5974 

- 0.0381 

- 0.5941 

0.5784 

- 0.0407 

- 0.5750 

0.5594 

- 0.0433 

- 0.5559 

0.5405 

- 0.0457 

- 0.5369 

0.5216 

- 0.0481 

- 0.5181 

0.5029 

- 0.0502 

- 0.4994 

0.4843 

- 0.0522 

- 0.4809 

0.4659 

- 0.0540 

- 0.4626 

0.4476 

- 0.0555 

- 0.4444 

0.4295 

- 0.0568 

- 0.4264 

0.4115 

- 0.0578 

- 0.4086 

0.3937 

- 0.0586 

- 0.3909 

0.3761 

- 0.0590 

- 0.3734 

0.3586 

- 0.0592 

- 0.3560 

0.3412 

- 0.0592 

- 0.3388 

0.3240 

- 0.0589 

- 0.3218 

0.3070 

- 0.0584 

- 0.3049 

0.2901 

- 0.0577 

- 0.2883 

0.2735 

- 0.0569 

- 0.2718 

0.2571 

- 0.0559 

- 0.2557 

0.2410 

- 0.0549 

- 0.2398 

0.2251 

- 0.0537 

- 0.2243 

0.2097 

- 0.0524 

- 0.2091 

0.1946 

- 0.0510 

- 0.1943 

0.1799 

- 0.0495 

- 0.1800 

0.1656 

- 0.0479 

- 0.1661 

0.1518 

- 0.0463 

- 0.1526 

0.1384 

- 0.0446 

- 0.1396 

0.1256 

- 0.0428 

- 0.1272 

0.1132 

- 0.0409 

- 0.1152 

0.1015 

- 0.0390 

- 0.1038 

0.0902 

- 0.0370 

- 0.0930 

0.0796 

- 0.0349 

- 0.0827 

0.0696 

- 0.0328 

- 0.0731 

0.0602 

- 0.0307 

- 0.0641 

0.0514 

- 0.0285 

- 0.0556 

0.0433 

- 0.0263 

- 0.0478 

0.0358 

- 0,0240 

- 0.0407 

0.0291 

- 0.0217 

- 0.0342 

0.0230 

- 0.0193 

- 0.0283 

0.0176 

- 0.0169 

- 0.0230 

0.0129 

- 0.0145 

- 0.0183 

0.0089 

- 0.0121 

- 0.0143 

0.0057 

- 0.0097 

- 0.0108 

0.0031 

- 0.0072 

- 0.0077 

0.0014 

- 0.0048 

- 0.0050 

0.0003 

- 0.0023 

- 0.0025 

0.0000 

0.0002 

0.0 

0.0004 

0.0027 

0.0027 

0.0016 

0.0051 

0.0058 

0.0035 

0.0076 

0.0094 

0.0061 

0.0100 

0.0135 

0.0095 

0.0124 

0.0183 

0.0136 

0.0148 

0.0237 

0.0184 

0.0171 

0.0297 

0.0240 

0.0194 

0.0364 

0.0303 

0.0217 

0.0437 

0.0373 

0.0240 

0.0518 

0.0450 

0.0262 

0.0604 

0.0534 

0.0285 

0.0697 

0.0624 

0.0307 

0.0796 

0.0721 

0.0328 

0.0902 

0.0824 

0.0350 

0.1013 

0.0933 

0.0371 

0.1129 

0.1047 

0.0391 

0.1251 

0.1168 

0.0411 

0.1379 

0.1294 

0.0431 

0.1512 

0.1425 

0.0449 

0.1649 

0.1562 

0.0467 

0.1791 

0.1703 

0.0485 

0.1938 

0.1849 

0.0501 

0.2089 

0.1999 

0.0517 

0.2245 

0.2154 

0.0531 

0.2404 

0.2312 

0.0545 

0.2567 

0.2474 

0.0558 

0.2733 

0.2640 

0.0570 

0.2902 

0.2809 

0.0581 

0.3074 

0.2981 

0.0591 

0.3249 

0.3156 

0.0600 

0.3427 

0.3333 

0.0607 

0.3606 

0.3512 

0.0614 

0.3787 

0.3694 

0.0619 

0.3970 

0.3876 

0.0623 

0.4155 

0.4061 

C . C 626 

0.4340 

0.4246 

0.0628 

0.4526 

0.4432 

0.0628 

0.4713 

0.4619 

0.0627 

0.4900 

0.4806 

0.0624 

0.5087 

0.4993 

0.0620 

0.5273 

0.5179 

0.0615 

0.5460 

0.5366 

0 . C 608 

0.5646 

0.5551 

0.0599 

0.5831 

0.5736 

0.0588 

0.6016 

0.5921 

0.0576 

0.6199 

0.6104 

0.0562 

0.6382 

0.6286 

0.0547 

0.6563 

0.6466 

0.0531 

0.6743 

0.6645 

0.0514 

0.6921 

C . 6 B 23 

0.0496 

0.7097 

0.6998 

0.0476 

0.7271 

0.7170 

0.0457 

0.7442 

0.7340 

0.0436 

0.7610 

0.7507 

0.0415 

0.7775 

0.7670 

0.0393 

0.7937 

0.7830 

0.0371 

0.8095 

0.7987 

0.0349 


- 0.6163 

- 0 . B 562 

1.0727 

- 0.6212 

- 0.8624 

1.0697 

- 0.6266 

- 0.8694 

1.0663 

- 0.6325 

- 0,8771 

1.0626 

- 0.6390 

- 0.8854 

1.0586 

- 0.6461 

- 0.8944 

1.0542 

- 0.6539 

- C .9044 

1.0492 

- C .6624 

- 0.9152 

1.0438 

- 0.6715 

- 0.9268 

1.0380 

- 0.6812 

- 0.9390 

1.0318 

- 0.6914 

- 0.9519 

1.0252 

- 0.7022 

- 0.9654 

1 . 0182 

- 0.7134 

- 0.9794 

1.0109 

- 0.7250 

- 0.9938 

1.0033 

- 0.7371 

- 1.0088 

0.9953 

- 0.7497 

- 1.0243 

0.9870 

- 0.7628 

- 1.0403 

0.9782 

- 0.7762 

- 1.0567 

0.9692 

- 0.7900 

- 1.0734 

0.9600 

- 0.8039 

- 1.0901 

0.9506 

- 0.8174 

- 1.1062 

0.9414 

- 0.8298 

- 1.1210 

0.9330 

- 0.8405 

- 1.1337 

0.9257 

- 0.8492 

- 1.1440 

0.9198 

- 0.8550 

- 1.1508 

0.9158 

- 0 . B 575 

- 1.1537 

0.9141 

- 0 . B 564 

- 1.1525 

0.9148 

- 0.8521 

- 1.1474 

0.9178 

- 0.8447 

- 1.1387 

0 . 922 B 

- 0.8350 

- 1.1272 

0.9294 

- 0.8238 

- 1.1140 

0.9370 

- 0.8121 

- 1.1000 

0.9450 

- 0.8006 

- 1.0862 

0.9528 

- 0.7896 

- 1.0729 

0.9602 

- 0.7791 

- 1.0602 

0.9673 

- 0.7689 

- 1.0478 

0.9741 

- 0.7590 

- 1.0356 

0.9808 

- 0.7489 

- 1.0233 

0,9875 

- 0.7386 

- 1.0106 

0.9943 

- 0.7281 

- 0.9976 

1.0013 

- 0.7172 

- 0.9841 

1.0084 

- 0.7059 

- 0.9700 

1.0158 

- 0.6940 

- 0.9552 

1.0235 

- 0.6813 

- 0.9392 

1.0317 

- 0.6676 

- 0.9218 

1,0405 

- 0.6531 

- 0.9034 

1.0497 

- 0.6376 

- 0.8835 

1.0595 

- 0.6205 

- 0.8616 

1.0700 

- 0.6018 

- 0.8374 

1.0815 

- 0.5810 

- 0.8104 

1.0940 

- 0.5578 

- 0.7798 

1.1077 

- 0.5310 

- 0.7445 

1.1230 

- 0.4997 

- 0.7027 

1.1404 

- 0.4625 

- 0.6526 

1.1600 

- 0.4179 

- 0.5920 

1.1822 

- 0.3636 

- 0.5172 

1.2068 

- 0.2960 

- 0.4229 

1.2334 

- 0.2092 

- 0.3001 

1.2604 

- 0.0955 

- 0.1374 

1.2823 

0.0545 

0.0785 

1.2862 

0.2439 

0.3494 

1.2506 

0.4599 

0.6492 

1.1614 

0.6723 

0.9278 

1.0375 

0.8504 

1.1454 

0.9189 

0.9830 

1.2967 

0.8282 

1.0774 

1.3987 

0.7643 

1.1408 

1.4644 

0.7224 

1.1796 

1.5035 

0.6972 

1.2027 

1.5263 

0.6825 

1.2140 

1.5374 

0.6753 

1.2184 

1.5417 

0.6725 

1.2181 

1.5414 

0.6727 

1.2160 

1.5394 

0.6740 

1.2137 

1.5371 

0.6755 

1.2116 

1.5351 

0,6768 

1.2097 

1.5332 

0,6780 

1.2082 

1.5317 

0,6790 

1.2072 

1.5308 

0.6796 

1.2068 

1.5304 

0.6798 

1.2069 

1.5305 

0.6798 

1.2075 

1.5311 

0.6794 

1.2086 

1.5322 

0.6787 

1.2102 

1.5337 

0.6777 

1.2120 

1.5355 

0.6765 

1.2141 

1.5375 

0.6752 

1.2163 

1.5397 

0.6738 

1.2189 

1.5422 

0.6722 

1.2216 

1.5448 

0.6705 

1.2246 

1.5478 

0.6686 

1.2276 

1.5507 

0.6667 

1.2308 

1.5537 

0.6647 

1.2342 

1.5571 

0.6626 

1.2380 

1.5607 

0.6602 

1.2420 

1.5646 

0.6577 

1.2462 

1.5686 

0.6550 
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RfiE 2822 CASE 9 SAMPLE CASE FOR LANGLET FINAL REPORT 
M=. 730 ALF= 2.706 CLC = 0.7977 NIT= 56 MESH=1C0X32 RSID=. 4230-09 
CL = 0.8031 CD =0.0166 COB =0.0166 CM =-.0875 CPTE=0. 1612 
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